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1 Einleitung 



Motivation dieser Arbeit ist ein Klassifikationsprojekt aus dem Gebiet der ganzzahli- 
gen Gitter. Martin Kneser beschrieb 1957 in dem Artikel Klassenzahlen definiter 
quadratischer Formen [Kne57] ein Exhaustionsverfahren für die Klassifikation voll- 
ständiger Geschlechter von Gittern. Diese Methode ist nicht unmittelbar praktikabel 
für ein automatisiertes Verfahren; sie läßt sich aber unter algorithmischen Aspekten 
modifizieren, um auch Geschlechter mit großer Klassenzahl klassifizieren zu können. 

Im Rahmen dieser Arbeit wurde daraus das Computerprogramm tn entwickelt, um 
die Geschlechter vieler Gitter zu bestimmen, die bis auf wenige Ausnahmen zuvor 
nicht vollständig bekannt waren. Bei der Konstruktion neuer Gitter mit Knesers Meth- 
ode behält man die Kontrolle über wichtige Gitterinvarianten. Damit unterscheidet 
sich tn von lokalen Optimierungs- oder Trainingsalgorithmen und hat eine gezielte 
Suche nach Gittern mit herausragenden Eigenschaften wie Extremalität, Wurzelsys- 
tem oder hoher Packungsdichte erlaubt. 

Die Algorithmik der Konstruktion ganzzahliger Gitter motivierte neue Problem- 
stellungen und erwies sich auch aus theoretischer Sicht als ein ergiebiges Gebiet. 
Das scheinbar banale Problem der Zerlegung eines Gitters in orthogonale Untergitter 
führte zur Entwicklung effizienter Algorithmen zum einen für das Zerlegungsproblem 
selbst, zum anderen für die Erzeugung eines Gitters aus einem großen Erzeugenden- 
system. Minkowskis zweiter Fundamentalsatz [Min96] aus der Geometrie der Zahlen 
führt zu Abschätzungen der Anzahl der benötigten arithmetischen Operationen und 
zu einer allgemeinen, quantitativen Aussage über die Größe minimaler Erzeugenden- 
systeme von Gittern. 

Der erste Teil der Arbeit leistet einen Beitrag zur konkreten Klassifikation ganz- 
zahliger Gitter. Die Fundierung und ein Teil der Ergebnisse wurden in der gemein- 
samen Arbeit Classification of integral lattices with large class number [SH98] mit 
Rudolf Scharlau publiziert. In diesem Kapitel wird auf die originären algorithmischen 
Aspekte der 2-Nachbarmethode fokussiert. Mit Hilfe von Knesers Nachbarverfahren 
gelang es H. Niemeier bereits 1973, das Geschlecht der geraden, unimodularen Gitter 
in Dimension 24 zu bestimmen [Nie73]. Er konnte dabei die starken Struktureigen- 
schaften der Gitter in diesem Geschlecht ausnutzen, um die Nachbargitter zu identi- 
fizieren und zu unterscheiden. In allgemeineren Situationen sind separierende Invari- 
anten aber nur mit großem algorithmischen Aufwand zu berechnen. In diesem Sinne 
äußern sich Conway und Sloane verhalten über den praktischen Nutzen von Knesers 
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Methode in nichttrivialen Fällen. 

„However there is a geometric method used by Witt and Kneser which 
(after the work of Niemeier) is effective roughly until the sum of the 
dimension and the (determinant)^/^ exceeds 24, beyond which point it 
seems that the forms are inherently unclassifiable." [CS92], S. 352 f. 

Diese Einschätzung erwies sich als zu pessimistisch, denn die Ergebnisse in [SH98] 
und in dieser Arbeit zeigen, daß Klassifikationen deutlich jenseits dieser Schranke 
durchführbar sind. 

Zwei verschiedene Gitter L und L' heißen benachbart, wenn sie ein gemeinsames 
Untergitter jeweils vom Index 2 enthalten. Bei geraden Gittern mit ungerader De- 
terminante läßt sich mit Knesers Nachbarmethode ein Nachbargitter L{v) aus L und 
einem Nachbarvektor v ^ L\2L mit Hilfe eines sogenannten Nachbarschrittes kon- 
struieren. Da wir nur an den Isomorphieklassen von Nachbargittem interessiert sind, 
läßt sich die Konstruktion aller Nachbargitter bis auf Isomorphie auf eine endliche 
Auswahl von Nachbarvektoren v zurückfuhren. Sie werden als geeignete Repräsen- 
tanten aus der endlichen Menge von Orbits der induzierten Automorphismengruppe 
von L auf L/2L ausgewählt. 

Das Computerprogramm tn basiert auf einer Implementation dieser Methode. Daneben 
enthält es zahlreiche Analysefunktionen zur Bestimmung von Invarianten der kon- 
struierten Gitter. Wegen der Komplexität der Berechnungen wurde das Design und 
die Realisierung von tn sehr stark an Performanceaspekten ausgerichtet und deshalb 
nicht in einem existierenden Computeralgebrasystem vorgenommen.^ 

Bei den enumerativen Klassifikationen konzentrierte sich das Interesse auf die £- 
modularen Gitter mit ^ = 3, 5, 7, 11 in Dimensionen bis zu 14 und die Bestimmung 
ihrer modularen und extremalen Gitter. Ein bemerkenswertes Ergebnis aus theoretis- 
cher Sicht ist ein Satz über das Geschlecht 7^. Die Klassifikation mit tn heferte den 
Beweis für die Nichtexistenz eines vermuteten extremalen Gitters. 

Das Programm tn kann zur Systematisierung der Ergebnisse eine Vielzahl von 
Invarianten eines Gitters berechnen. Eine bisher wenig beachtete Invariante ist die 
Fourier-Transformation der Längenfunktion eines Gitters auf L/2L. Erste Vermutun- 
gen, daß diese Fourier-Transformierte nur einen kleinen Träger besitzt, bewahrheit- 
eten sich nicht. Es scheinen jedoch im allgemeinen nur relativ wenig verschiedene 
Fourierkoeffizienten in der transformierten Funktion eine Rolle zu spielen. Mit tn 
lassen sich das Spektrum und die Fourier-Transformierte der Längenfunktion eines 
Gitters für eine feste Basis bestimmen. 

Die Nachbarschaftsmethode erweist sich auch in Fällen, in denen die Theorie keinen 
Erfolg garantiert, als erstaunlich robust. Bei Gittern mit gerader Determinante ist der 

'Seit 1999 existiert aucli eine Implementierang der Basisfunktionalität der 2-Nachbarmethode in dem 
Computeralgebrasystem MAGMA. 
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oben beschriebene Nachbarschritt nicht immer möglich. Dennoch läßt sich beispiel- 
sweise das (bereits bekannte) Geschlecht des Bames-Wall-Gitters [SV94] in der Di- 
mension 16 berechnen. Man kann beobachten, daß in höheren Dimensionen der Nach- 
barschaftsgraph, der gebildet wird aus den Isomorphieklassen von Gittern als Knoten 
und den Äquivalenzklassen von Nachbarvektoren als Kanten, für jedes Gitter so viele 
Kanten enthält, daß das Kollabieren einzelner Nachbarbildungen in der Regel nicht 
ergebnisrelevant ist. 

Das zweite Kapitel beginnt mit einer Diskussion des Zerlegbarkeitsproblems. Ein 
Gitter heißt zerlegbar, wenn es sich als direkte Summe nichttrivialer, orthogonaler 
Untergitter schreiben läßt. Die Zerlegung eines Gitters in unzerlegbare Gitter ist 
eindeutig bis auf die Reihenfolge der Summanden. Der konstruktive Beweis dieses 
Satzes stammt von Martin Kneser [Kne54]. Seine Konstruktion, eine Analogie zum 
Sieb von Erathostenes, dient primär dem Beweis der Behauptung, aber nicht der 
konkreten Durchführung. So erweist sie sich unter dem Blickwinkel der Berechen- 
barkeit als nicht ausreichend effizient. Wir zeigen in Satz 3.17, wie ein inkrementelles 
Konstruktionsverfahren statt der Benutzung des Siebverfahrens zu praktikablen Laufzeit- 
en führt, und geben eine Schranke für die Anzahl der algorithmischen Operationen des 
Gitterzerlegungsproblems an. 

Das Fehlen eines Basisergänzungssatzes und eines Steinitz'schen Basisaustauschsatzes 
in der Modulsituation verkompliziert die algorithmische Behandlung der Basisberech- 
nung eines Gitters außerordentlich. Das Hauptergebnis in diesem Kapitel ist ein 
Überdeckungssatz für Erzeugendensysteme von Gittern, der aus dem zweiten Fun- 
damentalsatz von Minkowski folgt. Es wird gezeigt, daß in typischen Situationen 
der praktischen Berechnung von Gittern ein mininales Erzeugendensystem nur aus 
wenig mehr Vektoren als dem Gitterrang besteht. Die quantitative Aussage liefert 
Satz 3.14: ist 5 C L ein minimales Erzeugendensystem eines Gitters L auf dem 
n-dimensionalen euklidischen Raum E mit Vektoren, die in der euklidischen Norm 
nicht länger als ß G M sind, dann ist 15 1 ^ n + log2(?^K5lET)")- ^^{e dieses 
Satzes wird gezeigt, daß das algorithmische Problem der Konstruktion von Gitter- 
basen aus großen Erzeugendensystemen durch das angegebene inkrementelle Ver- 
fahren effizient lösbar ist im Vergleich zu etablierten Methoden wie den Variationen 
des LLL- Algorithmus und der Berechnung der Hermite-Normalform. Tests mit ganz- 
zahligen Zufallsgittern unterstreichen die praktische Relevanz (siehe die Laufzeitver- 
gleiche auf Seite 45). 

Das dritte Kapitel behandelt die Berechnung von Gitterquantizern. Dabei handelt 
es sich um Approximationen, die aus der Informationstheorie motiviert werden. Ein 
Quantizer eines Gitters ordnet jedem Punkt im euklidischen Raum den nächstliegen- 
den Gitterpunkt zu. Er ist also im wesentlichen durch die Geometrie der einem Gitter- 
punkt nächstliegenden Punkte, der sogenannten Voronoi'zelle bestimmt. Die Qualität 
eines Gitterquantizers aus ingenieurwissenschaftlicher Sicht wird durch das dimen- 
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sionslose zweite normalisierte Trägheitsmoment der Voronoizelle (Definition 4.15) 
bewertet. Dies ist ein skaliertes Integral über die quadratische Vektorlängenfunktion 
in der Voronoizelle des Ursprungs. Die Abstandsfunktion wird algorithmisch durch 
das Verfahren von Fincke und Pohst auf Nebenklassen [FP85] ausgewertet.-^ Wir 
beschreiben die bekannten, grundlegenden Eigenschaften von Voronoizellen und die 
approximative Berechnung ihrer Gitterquantizer mit Monte-Carlo-Integration. Ihre 
Berechnung ist ebenfalls in tn implementiert. Alternativ lassen sich Gitterquantizer 
und Quantizer zu Kugelpackungskonstruktionen aus nichtlinearen Codes (Typ A Kon- 
struktion nach [CS92]) mit einem Programm für das Computeralgebrasystem MAG- 
MA berechnen. Gitterquantizer sind bis in die jüngste Vergangenheit nicht intensiv 
untersucht worden. In einer neueren Arbeit erzielten Agrell und Eriksson [AE98] eine 
Reihe von zum Teil kontraintuitiven Ergebnissen über Gitterquantizer. 

Für das Zustandekommen dieser Arbeit schulde ich vielen Personen Dank. An 
erster Stelle ist hier Herr Prof. Dr. Rudolf Scharlau zu nennen. Seine Betreuung 
gab dieser Arbeit den Anfang und die Richtung. Er war mir in all den Jahren ein 
aufmerksamer und geduldiger Gesprächspartner. Die Diskussionen mit meinem Kol- 
legen Frank Vallentin haben mir viele Anknüpfungspunkte an die Informatik eröffnet. 
Bernd Souvignier hat seine Implementationen der Berechnung von Erzeugern der Au- 
tomorphismengruppe eines Gitters und des Isometrietests zweier Gitter freundlicher- 
weise diesem Projekt zur Verfügung gestellt. Ohne die wunderbaren Entwicklungstools 
emacs, gcc und gbd aus dem GNU-Projekt von Richard Stallmans Free Software 
Foundation wäre aus diesem Projekt wahrscheinlich nicht ein so flexibles und leis- 
tungsfähiges Programm wie tn entstanden. Die Hochschulrechenzentren der Univer- 
sitäten Bielefeld, Dortmund und Duisburg haben mir Kapazitäten auf ihren Supercom- 
putern zur Verfügung gestellt. Der Besuch des AT&T Shannon Labs in New Jersey 
auf Einladung von Neil Sloane lenkte mein Interesse auf Gitterquantizer. Die Univer- 
sität Dortmund hat mir für die Anfertigung dieser Dissertation eine Assistentenstelle 
zur Verfügung gestellt. 

Ganz besonders danke ich meiner Frau Dr. Eva-Maria Kaffanke dafür, daß diese 
Arbeit nach langer Zeit doch noch ihre Vollendung gefunden hat. 



^In der Praxis von Telekommunikationsanwendungen wie bei Analog-Digital-Wandlern spielt neben 
der Güte des Trägheitsmomentes auch die Effizienz der Implementation der Approximation eine 
Rolle. 
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2.1 Einleitung 

In diesem Kapitel wird das Nachbarverfahren zur Klassifikation ganzzahliger Git- 
ter mit großer Klassenzahl beschrieben. Die zugrundeliegende Methode wurde von 
M. Kneser entwickelt [Kne57]. H.-V. Niemeier klassifizierte mit ihrer Hilfe 1973 die 
geraden unimodularen Gitter in Dimension 24 [Nie73]. 

Obwohl dieser Meilenstein in der Klassifikation der ganzzahligen quadratischen 
Formen die Bedeutung von Knesers Methode beweist, wurde der praktische Nutzen 
in „generischeren" Situationen als eher gering eingeschätzt (siehe [CS92], S. 352). In 
der gemeinsamen Arbeit mit R. Scharlau Classification of integral lattices with large 
class number [SH98] wurde gezeigt, daß eine algorithmisch optimierte Version des 
Nachbarverfahrens eine Vielzahl von Geschlechtem zugänglich macht, die mit bish- 
erigen Methoden nicht klassifizierbar waren. Besonders für die Suche nach Gittern 
mit bestimmten Eigenschaften kann ein Exhaustionsverfahren wie die Nachbarmeth- 
ode äußerst nützlich sein. Zu dieser Arbeit ist das Computerprogramm tn entstanden, 
dessen Funktionsweise im folgenden erläutert wird. In Abschnitt 2.2 werden die the- 
oretischen Grundlagen der Konstruktion von Gittern durch die 2-Nachbarmethode 
zusammengefaßt. In Abschnitt 2.5 wird die Funktion von tn genauer beschrieben 
und mit einigen markanten Resultaten in Abschnitt 2.7 illustriert. 

Knesers Verfahren garantiert unter gewissen Nebenbedingungen eine vollständige 
Klassifikation von (Spinor-)Geschlechtem. Die 2-Nachbarmethode läßt sich aber auch 
mit Erfolg bei Gittern mit gerader Determinante anwenden (Abschnitt 2.6), welche die 
Voraussetzungen des Satzes von Kneser nicht erfüllen. Ein interessantes Beispiel ist 
die Klassifikation des Geschlechtes des Barnes- Wall-Gitters (Abschnitt 2.7.2). 

2.2 Die Nachbarmethode 

Im folgenden seien alle Gitter ganzzahhg mit ungerader Determinante auf dem eu- 
klidischen n-dimensionalen Raum (V, (— , — )), sofern nicht etwas anderes vermerkt 
ist. 

Definition 2.1. Zwei Gitter L und L' werden Nachbarn genannt, falls 
L/{L n L') ^ Z/2Z ^ L'/{L n L'). 
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Aus dieser Definition folgt unmittelbar das 
Lemma 2.2. Zwei benachbarte Gitter L und L' besitzen dieselbe Determinante. 
Beweis. Nach der Determinanten-Index-Formel gilt 

detL = [L : L n L'f • det L n L' = [L' : L n L'f • detL n L' = detL'. □ 

Die effektive Konstruierbarkeit von Nachbarn klärt das 
Lemma 2.3. 

i) Für ein gerades Gitter L und einen Vektor v G L\2L ist das Gitter 

L(y) := Ly + ^Zv mit Ly := {x E L \ {x,v) G 2Z} 

ein Nachbar von L. Man nennt L{y) den Nachbarn von L bezüglich des Nach- 
barvektors V. 

ii) Wenn die Gitter L und L' Nachbarn sind, dann gilt L(2w) = L' für jedes 
weL'\L. 

Beweis. Zu i) Wir zeigen zunächst \L/Ly\ = 2. Die auf L := L/2L induzierte Bilin- 
earform (v.W) := {v,w) besitzt Diskriminante 1 und ist insbesondere regulär, weil 
die Diskriminante von (— , — ) ungerade ist. Wir nehmen an, daß L = Ly. Dann gilt 
{v, w) e 2Zfüi alle w e L. Somit ist {v, ItJ) = für alle w e L, also ist ?J = und 
damit v e 2L, was im Widerspruch zur Voraussetzung steht. 
Aus der Exaktheit der Sequenz 

0-^Ly-^L ^^M, Z/2Z 

folgt, daß |L/L^| =2. 

Dann gilt offenbar auch L/(LnL(f)) = 2. Da L gerade ist, folgt, daß v G L„,aber 
^ Ly, also ist L{v)/{L D L{v)) = Z/2Z und damit sind L und L{v) Nachbarn. 

Zuii). Es seiu) G L'\L für die Nachbarn L und L'. Dann folgt aus \L' /(LnL')\ = 
2, daß 2w e L\2L. Wie oben ergibt sich aus 

0^L2n,^L ^^(^^ Z/2Z 0, 

daß \L/L2w\ = 2. Weil L' ganzzahlig ist, gilt für jedes x G L', daß {x,2w) = 
2{x,w) G 2Z, also insbesondere L n L' C Da L/L n L' ^ Z/2Z ^ L/Ls^, 
muß \L2w/LnL'\ = 1 sein, also ist L2w = LHL'. Aus w G L'\Lund \L'/L2w\ = 2 
erhalten wir schüeßlich L' = L2w + \'^{'^w) = L{2w). □ 

Bemerkung 2.4. L{y) ist nicht notwendigerweise ganzzahlig. Dieser Sachverhalt wird 
im folgenden Lemma genauer untersucht. 
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2.3 Die Konstruktion des Nachibarscliaftsgraphen 

Lemma 2.5. Sei L ein gerades Gitter und v E L\2L. Dann ist 

i) L{v) ein ganzzahliges Gitter genau dann, wenn {v, v) G 4Z, 

ii) L{v) ist ein gerades Gitter genau dann, wenn {v, v) G 8Z. 
Beweis. 

i) Sei {v,v) G 4Z. Mit L ist auch ganzzahlig. Aus (L^,, i«) = C Z 
und {\v, \v) = \{v,v) C i4Z = Z, folgt, daß auch L{v) ganzzahlig ist. Die 
Umkehrung folgt unmittelbar aus der Ganzzahügkeit von L und {v^v) G Z <^=> 

ii) Sei {v, v) G 8Z. Mit L ist auch Ly gerade. Dann ist wegen (^■u, ^v) G | 8Z = 
2Z auch L(t;) gerade. Die Umkehrung folgt unmittelbar aus der Ganzzahügkeit 

von L und {v, v) e 2Z ^ {^v, Iv) G 8Z. 

□ 

Die beiden folgenden Lemmata erlauben es, die Konstruktion aller Nachbarn eines 
Gitters bis auf Isomorphie auf ein endliches Problem zurückzuführen. 

Lemma 2.6. Es seien w,w' G L \ 2L. Dann gilt L{w) = L{w'), falls w — w' £ 2L^. 

Beweis. Es sei z := G Lw. Dann folgt aus {x,w) = {x,w) — {x,2z) = 

{x, w — 2z) = {x, w') für alle x £ L, daß Lyj = Lyj'. Da 2L C = folgt 

L{w) = Lyj + ^Zw = Lyj + ^Z{w' + 2z) = L„> + ^Zw' = L{w'). 

□ 

Lemma 2.7. Sei L ein ganzzahliges, gerades Gitter und v G L\2L. Dann existiert 
ein y & L mit {v, y) ^ 2Z und für alle w & v & L/2L gilt entweder L{w) = L{v) 
oder L{w) = L{v + 2y). Weiter ist eines dieser Gitter gerade, das andere ungerade. 

Beweis. Da v ^ Q erfüllt sogar für jede gegebene Basis von L einer der Basisvek- 
toren die von y geforderte Eigenschaft. Ansonsten wäre die induzierte Form (— , — ) 
identisch und besäße damit nicht die Determinante 1. 

Falls {v, f ) = (mod 8), so ist L{v) ein gerades Gitter und L{v + 2y) wegen 

[v + 2y,v + 2y) = {v, v) + 4(t;, y) + 4(y, y) = A (mod 8) 

ein ungerades Gitter. Insbesondere sind L{v) und L{v + 2y) verschieden. Also gilt 
nach dem vorangehenden Lemma für alle w £ v wegen [2L : 2L^] = 2 entweder 
L{w) = L{y) oder L{w) = L{v + 2y). □ 



7 



2 Klassißkation ganzzahliger Gitter 



Damit braucht man lediglich 2" — 1 Repräsentanten der Klassen von L /2L = ¥2 
als Nachbarvektoren zu berücksichtigen, um alle Nachbargitter von L zu konstruieren. 
Dies reduziert die Komplexität der Auswahl der Nachbarvektoren auf ein endliches 
Problem. 

Dieses Argument beweist das 

KoroUar 2.8. Ein Gitter besitzt bis auf Isometrie höchstens 2" — 1 Nachbargitter. 

Das nächste Lemma zeigt, daß Nachbarvektoren, die in demselben Orbit unter der 
Automorphismengruppe des Gitters liegen, isometrische Nachbarn erzeugen. 

Lemma 2.9. Ist a G 0{V) und w e L\ 2L, dann ist L{wa) = L{w) a. 

Beweis. Für x e L gilt 

X & Lyjcr <J=> {x, wa) (x ,w) & 21, 4^ x G Lyj 4^ x G Ly, a. 

□ 

Bemerkung 2.10. Aufgrund des vorhergehenden Lemmas können wir uns bei der 
Auswahl von Nachbarvektoren auf diejenigen Gittervektoren beschränken, die ein- 
erseits in verschiedenen Orbits unter der Automorphismengruppe 0{L) liegen und 
die andererseits in verschiedenen Klassen modulo 2L hegen. Das kommutative Di- 
agramm in Abbildung 2.1 veranschaulicht diesen Sachverhalt. Alle enthaltenen Ab- 
bildungen sind surjektiv. Wegen der Kommutativität des Diagramms kann man die 
Berechnung aller Nachbargitter (bis auf Isometrie) von L vereinfachen, indem man 
zunächst zu einem Repräsentantensystem von L/2L = übergeht und darauf die 
induzierte Automorphismengruppe 0{L) in kanonischer Weise operieren läßt. Hier- 
bei bildet man Vektoren aus L nach F2 ab, vermöge Y17=i ' ' ' 
eine beliebige, aber feste Basis {fei, • • • , bn} von L. 



L 



¥0 



0{L) 



{vO{L) \ v€L} 



0{L) 



Abbildung 2.1: Die Orbits unter 0{L) 



Die Konstruktion aller Isomorphieklassen von Nachbargittern des Gitters L ist fol- 
gendermaßen möghch: 
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i) Wähle eine beliebige, von nun an fixierte Basis von L. 

ii) Identifiziere L/2L mit F2 und berechne alle Orbits von 0{L) auf F2 bezüglich 
der gewählten Basis. 

iü) Wähle aus jedem Orbit in unter der Operation von 0{L) einen Vektor v G 
F2 und suche ein beüebiges Urbild v e L auf. 

iv) Bilde den geraden Nachbarn von L bezüglich v, sofern er überhaupt existiert. 

Die Isomorphieklasse des Nachbargitters hängt nicht vom gewählten Repräsentanten 
ab, deshalb sprechen wir auch gelegentlich von einer Nachbarbildung „bezüglich der 
Klasse v" anstatt „bezüglich des Vektors v". 

Mit der iterierten Bildung von Nachbarn eines gewählten Gitters werden weitere 
Gitter konstruiert. Bei geraden Gittern besteht diese Menge aus einer Vereinigung 
von echten Spinorgeschlechtern, bei allen Beispielen und Resultaten in diesem Kapi- 
tel sogar aus einem ganzen Geschlecht (siehe [SH98], S. 742; [0'M71], §102 und 
[Kne56], Satz 2 bis Satz 5). Es bietet sich an, die Nachbarschaft von Gittern als Inzi- 
denzrelation in einem Graphen zu beschreiben. Dieser Graph ist ungerichtet, weil 
nach Lemma 2.3 ii) Nachbarschaft eine symmetrische Relation ist. 

Da wir im folgenden an einer Klassifikation dieser Geschlechter interessiert sind, 
definieren wir den Nachbarschaftsgraphen nicht auf der Menge aller Gitter, sondern 
nur innerhalb eines Geschlechts und nur bis auf Isomorphie der Gitter. 

Definition 2.11. Für ein Gitter L sei [L] die Klasse aller zu L isometrischen Gitter. 
Dann nennen wir den ungerichteten Graphen {V,E) mit V = {[L\ \ L ^ Q} und 

E = { {[Li], [-£'2]} I Li, L2 G Li und L2 sind Nachbargitter} 

den Nachbarschaftsgraphen von Q. 

2.4 Konstruktion einer Basis eines Nachbarn 

Sei {61, • ■ • , bn} eine Basis des geraden Gitters L und f G L \ 2L. In Algorithmus 1 
wird eine Basis des geraden Nachbarn bezüglich v konstruiert. Im Anschluß zeigen 
wir die Korrektheit des Verfahrens. 

In der Initialisierung im ersten Schritt wird v, wenn möglich, auf eine durch 8 
teilbare Norm gebracht. Ist {v,v) = (med 8), so bilden wir den Nachbarn L{v). 
Ist {v,v) = 4 (med 8), so wählen wir nach Lemma 2.7 ein bi mit {v,bi) ^ 2Z und 
bilden ^(^ + 26^). 
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Wir suchen den ersten ungeraden Koeffizienten vj. in v auf und multiplizieren v mit 
einem geeigneten Skalar x, so daß x ■ Vk = l (med 8) . In der Implementation wird 
X durch einen einfachen Table-Lookup bestimmt. 

3, falls Vk = 3 
X = 5, falls Vfc = 5 . 
7, falls Vk = 7 

Da alle skalaren Faktoren ungerade sind, gilt immer noch L{v) = L(xv) . Anschließend 
wird ein m ^ k bestimmt mit {v,hm) ^ 2Z. Angenommen, ein solches m existiere 
nicht. Sei M = {{hi}i^k)i- Dann liegt nach der Annahme v G M , also ergibt sich 
beim Übergang zum orthogonalen Komplement M C y-^. Aus Dimensionsgründen 
gilt M = U-L. Da {v, v) e 2Z, folgt v Ev-^ =Ji. Wegen ü^fc = T folgt v ^ M, was 
ein Widerspruch ist. 

Aus der Konstruktion folgt, daß b'^ G L^, für i / k und daß 5'^ G L{v), also 
gilt {b[,--- , b'n)z C L{v). Die Übergangsmatrix von • • • , 6„} zu {b[,--- , b'^} 
besitzt folgende Gestalt (hierbei nehmen wir ohne Beschränkung der Allgemeinheit 
an, daß k = 1 und m = n ist): 



1 


• • 


• ^\ 


* 1 


• • 


• 


* 


1 • 


• 




1 • 





* • 1 
\* * * • * * 2/ 



Die Determinante dieser Übergangsmatrix ist gleich 1 ; damit folgt aus dem Vergle- 
ich der Determinanten det{bi, • • • , b'^)z = det L = det L{v), daß {b'^, - ■ ■ , 6^)z = 
L{v). 



2.4.1 Ungerade Gitter und der duale Nachbarschaftsgraph 

Nach Lemma 2.7 gehört zu jedem geraden Nachbarn L' eines geraden Gitters L 
genau ein ungerader Nachbar L". Dieser Sachverhalt läßt sich auch von einem ele- 
mentaren algebraischen Standpunkt aus betrachten. Sei v e L \ 2L mit L{v) = L' 
und 2/ G L mit {v,y) ^ 2Z und {v,v) = (med 4), d.h., L{v + 2y) = L" . Wh" 
setzen M = L + ^Zv und erhalten mit M/L^ die elementarabelsche Gruppe der 
Ordnung 4 mit den drei nichttrivialen Elementen L/Ly, L' /Ly, L" /Ly. Mit L und v 
ist die Gruppe M/ Ly und damit das Paar aus dem ungeraden und aus dem geraden 
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Nachbarn von L eindeutig bestimmt. Dies liefert uns eine (wohlbekannte) Interpre- 
tation des Nachbarschaftsgraphen. Eine Kante im Graphen kennzeichnet eine Nach- 
barschaft von geraden Gittern. Dieser Kante ordnen wir das entsprechende ungerade 
Gitter zu. Der duale Graph ist also eine Teilmenge des Nachbarschaftsgraphen für 
ungerade Gitter. Bei ungeraden Gittern gibt es im allgemeinen keine Korrespondenz 
von geraden und ungeraden Nachbarn, wie man leicht an dem Fall sieht, daß v in dem 
orthogonalen Komplement eines ungeraden Teilgitters N c L liegt. Dann ist auch 
N C Ly, also sind beide Nachbarn von L bezügüch v ungerade. 

In Spezialfällen kann der Übergang zum dualen Graphen äußerst nützlich sein. 
R.E. Borcherds benutzt in [Bor92] den Nachbarschaftsgraphen der Niemeier-Gitter, 
um die ungeraden unimodularen Gitter in Dimension 24 zu klassifizieren. 

2.5 Das Programm tn 

Die konkrete Berechnung des Nachbarschaftsgraphen erweist sich als eine komplexe 
Aufgabe. Dies liegt an der Vielzahl der involvierten Teilprobleme und verwendeten 
Zahlsysteme. Im Programm selbst werden drei verschiedene Arithmetiken eingeset- 
zt. Der aufwendigste Teil in den Berechnungen von tn besteht in der Auswahl der 
geeigneten Vektoren zur Nachbarbildung. Das Implementationsziel von tn ist, alle 
Nachbarbildungen berechnen zu können, solange es überhaupt möglich ist, die Auto- 
morphismengruppen der beteiligten Gitter zu berechnen und die Gitter auf Isometrie 
zu prüfen. 

Das Programm tn behandelt nur ganzzahlige Gitter, so daß die meisten Berech- 
nungen in hardwarenaher 32-Bit Integerarithmetik erfolgen können. 

Einige Routinen verlangen die Berechnung der „kurzen Vektoren" eines Gitters; 
damit ist eine Auflistung oder Abzählung aller Gittervektoren bis zu einer bestimmten 
Länge, meist der Länge des größten Basisvektors, gemeint. Dazu benutzen wir den 
Algorithmus von U. Fincke und M. Pohst [FP85] mit den Ergänzungen von H. Cohen 
(Algorithm 2.7.7 in [Coh96], S. 105) in einer Implementation mit Doublearithmetik. 

Die Berechnungen zur Auswahl der Vektoren zur Nachbarbildung werden, wie in 
Bemerkung 2.10 beschrieben, hauptsächlich im Vektorraum L/2L = ¥2 vorgenom- 
men. Für diesen Zahlbereich wurde eine eigene Implementation in Binärarithmetik 
implementiert, die in Abschnitt 2.5.2 vorgestellt wird. 

2.5.1 Die Berechnung der Automorphismengruppe eines Gitters 

Für die Berechnung von Automorphismengruppen und Isometrieen von Gitteren wer- 
den in tn die Programme autom und isom verwendet, die uns Bernd Souvignier 
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freundlicherweise zur Verfügung stellte. Sie sind Implementation der von B. Souvi- 
gnier und W. Plesken in [PS97] beschriebenen Verbesserungen der generischen Sims- 
StabiHsatorkettenmethode in Gittern. 

Die Routinen autom und isom sind vom Rechenzeitbedarf her die aufwendigsten 
Programmteile von tn. 



2.5.2 Rechnen in L/2L 

Da L = L/2L = F2 , hegt es nahe, die Rechnungen in L in hardwarenaher Binärarith- 
metik zu implementieren. Die Berechnung aller Orbits unter 0{L) erfolgt durch An- 
wendung aller von autom berechneten und anschließend induzierten Erzeuger. Aus 
diesem Grund ist die Effizienz dieser Routinen im Hinblick auf Geschwindigkeit und 
Speicherplatzbedarf in der Praxis ein kritischer Erfolgsfaktor des Programms tn. Die 
Implementation erfolgte in der Programmiersprache „C" [KR88]. 

Vektoren aus F2 werden mit dem Datentyp unsigned int identifiziert, der in 
den meisten Prozessorarchitekturen in einem 32-Bit Wort implementiert ist. In Di- 
mensionen jenseits von 32 sind aus mehreren Gründen keine kompletten Klassifikati- 
onen zu erwarten. Auf der einen Seite dauern Automorphismengruppenberechnungen 
und Isometrietests bereits sehr lange, auf der anderen Seite ist die Klassenzahl bei 
Geschlechtern in diesen Dimensionen sehr groß. Beispielsweise weiß man durch An- 
wendung der Maßformel (s.u.), daß es mehr als achtzig Milhonen unimodulare gerade 
Gitter in Dimension 32 gibt. 

Bei der Berechnung der Orbits von 0{L) auf L verwenden wir zwei verschiedene 
Kodierungen für Vektoren in F", je nach Optimierungsziel, Rechenzeitbedarf und 
Speicherplatzbedarf. Bei der Matrix- Vektor-Multiplikation wird ein Vektor von L 
als Bitfolge in einem Objekt vom Typ unsigned int kodiert. Eine Matrix besteht 
aus einem Array von n solcher Objekte. Für die Markierung der Orbits von 0{L) 
verwenden wir im Vektorraum L ein Kompressionsverfahren, bei dem jeder Vektor in 
L nur durch ein einziges Bit repräsentieren wird, und die Kodierung und Dekodierung 
sehr effizient erfolgt. Die folgenden beiden Beispiele illustrieren die benutzten Tech- 
niken bei der Matrix- Vektor Multiplikation in L (2.5.2.1) und der Bestimmung eines 
kanonischen Repräsentanten eines Orbits von 0{L). Der Quelltext von tn ist hier für 
den Zweck der Darstellung aufbereitet worden. Fehlerbehandlung und der Umgang 
mit Spezialfällen wie dim L ^ 3 führen zu etwas komphzierterem Code. 



'Bernd Souvignier implementierte später neben isom und autom die 2-Nachbarmethode in dem 
Computeralgebrasystem MAGMA. 
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2.5.2.1 Matrix-Vektor-Multiplikation in L 

Als Beispiel für die Effizienz der Arithmetik in L geben wir hier die verwendete Rou- 
tine für eine Matrix- Vektor-Multiplikation an. 



1 /* Matrix — Vektor Multiplikation in L/2L */ 
2 

3 typedef unsigned int vector_2 ; 
4 

5 /* A definition for matrices over GF(2). */ 

6 typedef 

7 struct struct_matrix_2 

8 ( 

9 int dim; /* dimension */ 

10 vector_2 *v; /* rows of the matrix */ 

11 } 

12 *matrix_2; 
13 

14 vold 

15 matrix_vector_2 (vector_2 *imv , const matrix_2 m, 

16 const vectür_2 v) 

17 /* matrix_veclor_2 compules ihe image of */ 

18 /* the column veclor v und er the */ 

19 /* Operation of the matrix m and puls the */ 

20 /* result in imv . */ 

21 ( 

22 int i , j ; 

23 vector_2 tmp; 
24 

25 for (i = 0, *imv = 0; i < m— >dim; i++){ 

26 /* Compute (m—>v[ij,v) and set result 

27 in f—!h hft of im */ 

28 tmp = m— >v [ i J & v ; 

29 for (j = _WORDSIZE / 2; j >= 1 ; j /= 2) 

30 tmp ^= tmp » j ; 

31 *imv 1= (tmp & 1) « ((m— xlim — 1) — i); 

32 ) 

33 ) 



Listing 2.1: Effiziente Matrix- Vektor Multiplikation in L 

Nachdem der Bildvektor imv mit dem Nullvektor initialisiert ist, wird in der f or- 
Schleife jeweils die i-te Komponente von imv ermittelt. Hierzu berechnet man das 

Skalarprodukt der i-ten Zeile m->v [i] der Matrix m mit dem Vektor v und weist 
es dem i-ten Bit in im zu. Man setzt j auf die halbe Anzahl der Bits in einem 
vector_2. 

Die Zeile 30 wird in der Schleife durchlaufen, um durch mehrfache Faltung mit 
einer XOR-Operation das Ergebnis im 0-ten Bit von tmp zu erhalten. Die oberen Bits 
spielen für das Ergebnis dann keine Rolle mehr. Das 0-te Bit wird in Zeile 3 1 an die 
(n — i)-te Stelle in im kopiert. Um bereits gesetzte Bits nicht zu löschen, erfolgt dies 
mit einer logischen OR-Operation auf im . 

2.5.2.2 Orbits unter Gruppenoperation 

Die allgemeine Standardmethode, einen Orbit unter einer Gruppenoperation zu bes- 
timmen, ist der elegante union-find Algorithmus mit Pfadkompression [TvL84]. 
In unserer Situation können wir zusätzlich ausnutzen, daß die Orbits aller Vektoren in 
L bestimmt werden müssen. Dies erlaubt ein effizienteres Vorgehen sowohl in Bezug 
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auf Speicherplatz als auch auf Rechenzeit. Dabei identifizieren wir L mit den Koordi- 
natenvektoren, die durch eine feste Basis von L gegeben sind und kodieren jeden Ko- 
ordinatenvektor auf ein bestimmtes Bit. Diese Bits bilden eine Folge in einem zusam- 
menhängenden Speicherbereich. Dadurch wird auf L eine totale Ordnung induziert, 
die bei unserer Kodierung mit der üblichen lexikographischen Ordnung zusammen- 
fällt. Dies erlaubt eine sehr schnelle Kodierung und Dekodierung von Vektoren in L. 
In tn sind für eine Liste von Vektoren in L die folgenden Methoden für Listenopera- 
tionen implementiert. 

insert_vector (v) Füge den Vektor v in die Liste ein. 
delete_vector (v) Lösche den Vektor v aus der Liste. 

lookup_vector (v) Prüfe, ob der Vektor sich bereits in der Liste befindet und 
liefere dementsprechend true oder f alse zurück. 

f ind_min Liefere den bezüghch der lexikographischen Ordnung kleinsten Vektor 

in der Liste zurück. 

Diese Operationen sind aus Effizienzgründen in tn zum Teil inline implementiert, d.h. 
nicht durch explizite Funktionsaufrufe. Der Klarheit halber formulieren wir beispiel- 
haft zwei dieser Operationen als Funktionsaufrufe. 

Die Vektoren werden in einem Array von unsigned char, einem 8-Bit Daten- 
typ kodiert. Dabei geben die führenden n — 3 Bits die Position des Bytes im Array 
an, in welchem der Vektor kodiert wird, und die letzten drei Bits werden durch die 
Position in dem entsprechenden Byte kodiert. 

Als erstes Beispiel beschreiben wir die Operation insert_vector ( v) . Um das 
den Vektor v beschreibende Bit aufzufinden, berechnen wir zunächst den Offset im 
Array. Die führenden n — 3 Bit von v, interpretiert als Binärzahl, geben die Position 
des entsprechenden Bytes im Array an (im Source Zeile 14 1 [ v»3 ] ). In diesem Byte 
setzen wir das Bit, dessen Position durch das Bitmuster letzten 3 Bit von v in der 
üblichen Binararithmetik beschrieben wird. Hierzu maskieren wir alle Bits von v bis 
auf die letzten 3 aus ( ( v & 7 ) ) und verwenden das Ergebnis, um das erste Bit um 
die entsprechende Anzahl zu verschieben (im Source: 1 << (v & 7)). 

1 /* insert vector v in list l */ 

2 

3 typedef unsigned int vector_2 ; 

4 typedef *unsigned cliar vlist; 
5 

6 /* Insert v in list I . Idea : Lookup ihe hyle 

7 in vlist l whose offset is binary coäed by 

8 the leading n—3 bits in v. 

9 Now mark in this byte the bit whose number 
10 is hinary coded by the last 3 bits. */ 

11 

12 void 

13 insert_vector (vlist 1, vector_2 v){ 

14 l[v>>3] 1= (1 « (v & 7))\ label { insert_vector_fast_decoding } 
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Listing 2.2: Operationen in L: insert_vector (v) 

B ei der folgenden Routine f i n d_m i n ( L ) , die einen Vektor aus der Liste L zurück- 
liefert, müssen wir die Fehlersituation der leeren Liste abfangen. Der Nullvektor spielt 
in L keine Rolle, weil er nicht zur Nachbarbildung verwendet werden kann. Wir be- 
nutzen deshalb sein Kodifikat „0" als Flag, um der aufrufenden Funktion mitzuteilen, 

daß die Liste leer ist. 



1 


int dim ; /* // 


e dimension of \overline {Lj *■/ 


2 


typedef unslgned int vector_2 ; 


3 


typedef *unsigned cliar vlist; 


4 


/* Lookup the 


"smallest" vector in l 


5 
6 


and return 


it . */ 


7 


vector_2 




8 


find_min (vlist 1){ 


9 


int i , t ; 




10 


/* Lookup the 


firsl non—empty byte 


11 


in the array L */ 


12 


for (i = 0; 


i < (1 « (dim - 3)); i++) 


13 


if (l[i] 


= 0) 


14 


brealt ; 




15 






16 


for (t = 0; 


t < 8; t++) 


17 


if (((l[i] 


» t) & 1) == 1) 


18 


brealt ; 




19 


if (t == 8) 




20 


return 0; 


/* There is no umarked vector. */ 


21 


eise 




22 


return (8 


* i + t ); 


23 


) 





Listing 2.3: Operationen in L: f ind_min 



Um in der Liste den kleinsten Vektor aufzufinden, d.h., in dem Array das erste 
gesetzte Bit zu finden, suchen wir die Liste linear durch. Diese Implementation ist 
für unsere Zwecke in der Laufzeit hinreichend effizient, denn wir können jeweils 
8 Bits zugleich prüfen, in dem wir jeweils ein ganzes Byte betrachten und testen, 
ob es gleich ist.^ Beispielsweise muß in der Dimension 20 lediglich ein zusam- 
menhängender Speicherbereich von höchstens 128 kByte durchlaufen werden. Der 
damit verbundene Aufwand ist irrelevant gegenüber den anderen Operationen in tn 
wie beispielsweise der Basisreduktion. Dieses Durchlaufen des Arrays geschieht in 
der folgenden Schleife. 

for (i = 0; i < (1 << (dim -3)); i++) 
if (l[i] != 0) 
break; 

Nun enthält die Variable i den richtigen Offset und wir suchen in dem i-ten Byte des 

Arrays das erste gesetzte Bit. 

^Die naheliegende Verbesserung, gleich 32 oder mehr Bits durch einen Gast auf einen größeren Daten- 
typ wie Integer zu überprüfen, führt nicht zu wesentlichen Geschwindigkeitssteigerungen und macht 
die Portierung auf andere Betriebssysteme komplizierter. 
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for (t = 0; t < 8; t++) 
if ( (1 [i] » t) & 1) 
break; 

Als Ergebnis geben wir den Vektor zurück, dessen Bitmuster der Binärzahl 8i + t 
entspricht. 

2.5.3 Die diskrete Fourier-Transformierte der Längenfunktion 

2.5.3.1 Die diskrete Fourier-Transformation auf endiicfien abeiscfien 
Gruppen 

Das Programm t n enthält eine Vielzahl von Analysefunktionen, um Invarianten eines 
Gitters zu bestimmen. Darunter befindet sich die bisher wenig untersuchte Funktion 
•^L/2L (^) ! die Fourier-Transformierte der Längenfunktion £. Diese Zerlegung der Län- 
genfunktion, welche einer Nebenklasse in L/2L die Norm eines kürzesten in ihr en- 
thaltenen Vektors zuordnet, in eine Linearkombination aus irreduziblen Charakteren 
von L/2L erlaubt bei vielen Gittern eine recht einfache Beschreibung dieser Funk- 
tion. Die Zerlegung selbst ist abhängig von der Wahl der zugrundeliegenden Basis in 
L beziehungsweise in L/2L; die Vielfachheiten der Koeffizienten in der Linearkom- 
bination, das sogenannte Spektrum, sind hingegen unabhängig von der Wahl einer 
Basis. Wir beschreiben zunächst die Fourier-Transformation und ergänzen diese mit 
einigen Beispielen. 

Sei G eine endliche abelsche Gruppe, die Menge der irreduziblen Charaktere 
von G und seien /, /' komplexwertige Funktionen auf G. Die irreduziblen Charaktere 
G^ bilden eine Orthonormalbasis bezüglich des Skalarproduktes (/, /') = X^geG fi9)f'{9) 
im Raum der Klassenfunktionen G* auf G. Deshalb kann man / folgendermaßen als 
Linearkombination der irreduziblen Charaktere schreiben, die eine Orthonormalbasis 
im Raum der Klassenfunktionen bilden: 

X6GV 

Man nennt /(x) := {x, f) = py Egec den Fourierkoeffizienten von x 

bezüglich /. 

Dadurch wird der Operator JT der Fourier-Transformation auf der Gruppe G definiert 

J^g: G* (G*)* 

f f ■x^<xJ > 

Wenn man die Charaktertafel von G als Matrix (x(5'))xeGv,geG und / als Spal- 
tenvektor {f{g))g^G schreibt, dann berechnet sich (/(x))x ^^f einen Faktor als 
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einfache Matrizenmultiplikation durch 

(/(X))x = 1^ • ix{9))x,9 U{9))g ■ 

2.5.3.2 Die diskrete Fourier-Transformierte der Längenfunl<tion eines 
Gitters 

Für ein Gitter sei l : L/2L — ^ R mit i{v) = min^^y^w, w) als die Längenfunktion 
von L definiert. 

2.5.3.2.1 Der Operator /2L Die Gruppe L/2L ist eine elementarabelsche Gruppe 
der Ordnung 2" mit n = rankL. Im ersten Schritt verifiziert man, daß sich ihre 
Charaktertafel in Form einer Hadamard-Matrix schreiben läßt. Diese Matrix ist bis 

auf einen skalaren Faktor gleich der Operatormatrix von Tj^j^L- Zunächst benötigt 
man eine Notation für die irreduziblen Charaktere von L/2L = Die Charak- 
tere dieser Gruppe lassen sich als mehrfache Tensorprodukte der Charaktere von F2 
schreiben. 

Die Charaktertafel von F2 enthält genau zwei Charaktere, den Hauptcharakter 
und den Signumcharakter x' . 








1 




1 


1 


X 


1 


-1 



Tabelle 2.1: Die Charaktertafel von F2 

Die irreduziblen Charaktere von F2 ergeben sich als n-fache Tensorprodukte der 
Charaktere x"*" und x~- Zur Vereinfachung der Notation wird ein irreduzibler Charak- 
ter von F" nur mit einer Folge von -|- und — indiziert und nicht als voll ausgeschriebenes 
Tensorprodukt, beispielsweise x^~^^ ■= <^ X~ ® <S) x^- Damit ergeben sich 
für F2n und F| und F^ die Charaktertafeln wie sie in den Tabellen 2. 1-2.3 beschrieben 
sind. 

Das Programm tn berechnet die Fourier-Transformierte der Längenfunktion £. 
Bei der Berechnung einer Vielzahl von Gittern stellt sich heraus, daß sich die Koef- 
fizienten in der Zerlegung in irreduzible Charaktere auf nur wenige Werte bei großen 
Vielfachheiten konzentrieren. 

Die folgenden Beispiele sollen diese Beobachtung illustrieren und zu neuen Fragestel- 
lungen Anlaß geben. 
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00 


Ol 


10 


11 




1 


1 


1 


1 




1 


-1 


1 


-1 




1 


1 


-1 


-1 


X 


1 


-1 


-1 


1 



Tabelle 2.2: Die Charaktertafel von F| 





000 


010 


100 


110 


001 


011 


101 


III 




1 


1 


1 


1 


1 


1 


1 


1 




1 


-1 


1 


-1 


1 


-1 


1 


-1 




1 


1 


-1 


-1 


1 


1 


-1 


-1 


x+" 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


X-++ 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


X-+- 


1 


-1 


1 


-1 


-1 


1 


-1 


1 


X-+ 


1 


1 


-1 


-1 


-1 


-1 


1 


1 


X 


1 


-1 


-1 


1 


-1 


1 


1 


-1 



Tabelle 2.3: Die Charaktertafel von F| 



2.5.3.3 Beispiele 

i) Sei L = (Z2, .)) mit G{q) = f J ^ ) ^ann ist 



(%)). = ^^-(x(5))x,P (%)).= 
LI 1 1 \ /0\ / 6 \ 





/l 1 1 


1\ 






/6\ 


1 


1 -1 1 


-1 


2 


1 


-4 


4 


1 1 -1 


-1 


1 


~ 4 


-2 




\l -1 -1 


1 / 


\^) 







Also läßt sich i in die folgende Linearkombination irreduzibler Charaktere zer- 
legen 

^ = ^(6X++ - 4x+- - 2x-+ + Ox~). 
ii) Sei A2 = .)) mit G{q) = ( ~o)- zerfällt ^ analog des 
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vorhergehenden Beispiels in die Linearkombination 



2x-+ - 2x^ 



2X" 



üi) Für größere n zeigt sich bei der Zerlegung von £ für die Gitter An eine Verteilung 
der Koeffizienten auf nur wenige Werte. 

a) Koeffizientenverteilung für A12 



Koeffizient 


Vielfachheit 


26624 


1 


-2048 


13 





4082 



b) Koeffizientenverteilung für ^13 



Koeffizient 


Vielfachheit 


57344 


1 


-4096 


14 





8177 



c) Koeffizientenverteilung für Au 



Koeffizient 


Vielfachheit 


122880 


1 


-8192 


15 





16358 



2.6 2-Nachbarn bei gerader Determinante 

Bisher haben wir die 2-Nachbarbildung nur bei geraden Gittern mit ungerader Deter- 
minante untersucht. Den Grund für diese Einschränkung finden wir im Beweis von 
Lemma 2.3. Ly muß bei der Nachbarbildung ein echtes Untergitter von L sein; dies 
ist im allgemeinen aber nicht gegeben, wie wir im Beispiel L = 2M für ein ganz- 
zahliges Gitter M sehen. Dann gilt für jeden Vektor v £ L, daß L = ist. Für viele 
Vektoren v ist L ^ aber auch bei gerader Determinante erfüllt und damit ist die 
Konstruktion eines Nachbarn ebenfalls möghch. 

Für die Konstruktion des 2-Nachbarschaftsgraphen eines Gitters L mit gerader De- 
terminante, das aus den Vektoren 61 , • • • ,bn erzeugt wird, geht man im wesentlichen 
genau so vor wie im ungeraden Fall. Man berechnet die Menge aller Orbits von 0{L) 
auf L und wählt aus jedem Orbit einen Repräsentanten v aus. Dann sucht man einen 
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Basisvektor bi von L mit {v,bi) = 1 und v ^ bi. Diese Eigenschaften sind unab- 
hängig von der Auswahl des Repräsentanten v. Wie in Lemma 2.3 i) gezeigt, ist bei 
ungerader Determinante die Existenz eines solchen 6j gesichert, im geraden Fall je- 
doch nicht. Findet man keinen Basisvektor 6j mit diesen Eigenschaften, dann verwirft 
man diesen Orbit und geht zum nächsten über. 

Dieses Vorgehen garantiert natürlich nicht eine erfolgreiche Klassifikation eines 
ganzen Geschlechtes, führt in der Regel aber in größeren Dimensionen (etwa ab Di- 
mension 5 bis 8) zum Ziel. In fast allen geraden Gittern L höherer Dimension finden 
sich genügend Orbits von Vektoren v ^ L mit L„ ^ L, welche die Existenz von 
2-Nachbarn oft sicherstellen. Wir verfolgen hier mit t n einen pragmatischen Ansatz, 
tn berechnet den Nachbarschaftsgraphen, und eine erfolgreiche Verifikation mit der 
Maßformel sichert die Vollständigkeit des Ergebnisses (s. Abschnitt 2.8). Wie oben 
bereits beschrieben, sind jedoch Beispiele von Gittern ohne 2-Nachbam auch in hohen 
Dimensionen konstruierbar. 

Als prominentes Beispiel für eine erfolgreiche Klassifikation in einem Geschlecht 
mit gerader Determinante führen wir in Abschnitt 2.7.2 die Klassifikation des Geschlecht- 
es des Barnes- Wall-Gitters an. 

2.7 Ergebnisse 

Die wichtigsten Ergebnisse, die mit Hilfe von t n erzielt wurden, sind in dem Artikel 
[SH98] dokumentiert. Wir geben hier eine kurze Zusammenfassung und einige neuere 
Resultate an. 

2.7.1 £-modulare Gitter 

Ein n-dimensionales ganzzahliges Gitter L heißt modular vom Level £, wenn es ein 
£ G N gibt mit L ^ \^,L* [Que95]. Es folgt unmittelbar, daß dct L = insbeson- 
dere ist n gerade, wenn £ quadratfrei ist. Viele der aus der Literatur bekannten Gitter 
sind modular, so z.B. das 12-dimensionale Coxeter-Todd-Gitter, das 16-dimensionale 
Bames-Wall-Gitter und die geraden unimodularen Gitter Eg,, das Gosset-Gitter, und 
A24, das Leech-Gitter [CS92]. Aus der Theorie der Modulformen ergibt sich für mod- 
ulare Gitter eine obere Schranke für das Minimum. Gitter, die dieses theoretische 
Maximum annehmen, heißen extremal. Für eine Einordnung und Diskussion dieses 
Themas verweisen wir auf den Übersichtsartikel [SSP99]. Die Existenz und die Ein- 
deutigkeit modularer Gitter ist eine zentrale Frage in dieser Theorie. Einige der ex- 
tremalen Gitter realisieren die besten bekannten Packungsdichten oder würden sie im 
Falle der Existenz übertreffen. 
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Für eine Berechnung des Nachbarschaftsgraphen mit tn benötigen wir ein Startgit- 
ter aus dem zu bestimmenden Geschlecht. Außer bei £ = 5 haben wir bei den unten 
aufgeführten Fällen (hier gilt jeweils 1 = 3 (mod 4)) mit dem direkten Produkt aus 

fi/2 l\ 

n/2 Gittern in der Dimension 2 mit Grammatrix ( ^ ^ ) begonnen. 

Wir zitieren hier eine Liste mit den aus der Theorie der Modulformen berechneten 
möglichen maximalen Minima. 



n 


4 


6 


8 


10 


12 


14 


16 


18 


£ = 


3 : 


min 


2 


2 


2 


2 


4 


4 


4 


4 


e = 


5 : 


min 


2 




4 




4 




6 




e = 


7 : 


min 


2 


4 


4 


4 


6 


6 


6 


8 


e = 


11 : 


min 


4 


4 


6 


6 


8 


8 







Tabelle 2.4: Maximal mögliche Minima von modularen Gittern (aus [SH98]) 



2.7.1.1 DerLevel^ = 3 

Die Klassenzahlen für die Dimensionen 2, 4, 6, 8, 10, 12 sind jeweils 1, 1, 1, 2, 3, 10; 
alle diese Gitter sind modular. Bis auf das Coxeter-Todd-Gitter in Dimension 12 
mit Minimum 4 sind alle anderen Gitter reflektiv, d.h., sie enthalten ein Wurzelsys- 
tem von vollem Rang ([SB96]). Eine Methode, das Geschlecht des Coxeter-Todd- 
Gitters über die Konstruktion der möglichen Wurzelsysteme zu bestimmen, und die 
daraus folgende Klassifikation wurden bereits in [SV95, SVOO] ohne Detailrechnun- 
gen angegeben. 

Mit dem Programm tn läßt sich aufgrund des Konstruktionsverfahrens nicht nur 
eine Auflistung aller Gitter eines Geschlechts berechnen, sondern auch der vollständi- 
ge Nachbarschaftsgraph beschreiben. Die folgende Tabelle enthält die Inzidenzmatrix 
des Nachbarschaftsgraphen des Geschlechtes des Coxeter-Todd-Gitters, bei welchem 
die Kanten mit der Anzahl der verschiedenen Orbits gewichtet wurden, die unter der 
jeweiligen Operation von 0{L) auf L/2L zu einer Nachbarschaftsbildung führen. 

In der Dimension 14 gibt es 29 Gitter in diesem Geschlecht, darunter ein extremales 
mit Minimum 4. Seine Automorphismengruppe von der Ordnung 2^ • 3^ • 7 • 13 ist 
(die modulo der Spiegelung durch den Ursprung einfache Gruppe) 02(3) • 2 ([SH98], 
Proposition 3.2). Eine explizite Konstruktion dieses Gitters wird im ATLAS [CCN+85], 
S. 60 beschrieben. 
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Kl K2 Ks K^ K5 Ke K^ Ks Kg K^q 



Kl 
K2 
Ks 
Ka 
K5 
Ke 
Kr 
Ks 
Kg 

KlQ 



1 1 



15 1111 

14 1 1 
117 1 III 
1 15 11 1 



1 III 



III 61 1 

I 14 1 

II 14 1 
1 1 2 



Tabelle 2.5: Der Nachbarschaftsgraph des Geschlechtes des Coxeter-Todd-Gitters mit 
den Bezeichnungen aus [SV95] 



2.7.1.2 DerLevel£ = 5 

Solche Gitter existieren nur in durch 4 teilbaren Dimensionen. In der Dimension 
12 existiert ein interessantes Geschlecht mit Determinante 5^, welches aus 48 Git- 
tern besteht, darunter 40 modularen und 4 extremalen. Die Klassifikation dieses 
Geschlechtes wurde unabhängig von G. Nebe (unveröffentlicht) erzielt. 



2.7.1.3 DerLevel£ = 7 

Das Hauptergebnis des Klassifikationsprojektes in [SH98] ist der Beweis der Nichtex- 
istenz eines vermuteten Gitters. In der Dimension 12 existiert kein extremales Gitter 
mit der Determinante 7^ . Dies ist das kleinste bekannte Beispiel für die Nichtexistenz 
eines möglichen extremalen Gitters. 



2.7.1.4 Der Level 11 

In den Dimensionen 4, 6, 8 existieren jeweils 3, 5, 31 Gitter, alle sind modular und in 
jedem Geschlecht gibt es ein eindeutiges extremales Gitter. 

In der Dimension 10 gibt es neben dem Craig-Gitter ein weiteres extremales (Beobach- 
tung von H.-G. Quebbemann, [Que95]) alle weiteren der insgesamt 297 Gitter be- 
sitzen ein kleineres Minimum. Siehe Abschnitt 2.7.3.1 für den 12-dimensionalen Fall. 
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2.7.2 Das Geschlecht des Barnes-Wall-G Itters Aie 



Das Gitter Aie wurde 1959 von Barnes und Wall konstruiert [BW59] und ist mit Min- 
imum 4 und Determinante 256 das dichteste bekannte Gitter in Dimension 16. Die 
Bezeichnung Aiq ist durch seine Konstruktion als geschichtetes Gitter motiviert (siehe 
[CS92], Chapt. 6). Das Geschlecht des Barnes-Wall-Gitters wurde von R. Scharlau 
und B.B. Venkov [SV94] klassifiziert. Analog zur Klassifikation der geraden, uni- 
modularen Gitter in Dimension 24 von B.B. Venkov in [Ven92] werden die Gitter im 
Geschlecht des Barnes-Wall-Gitters über ihr Wurzelsystem bestimmt. 

Obwohl es sich hier um ein Geschlecht von gerader Determinante handelt, set- 
zen wir tn mit den in Abschnitt 2.6 beschriebenen Modifikationen ein. Ausgehend 
von einer Erzeugermatrix für Aiq untersuchen wir zunächst beispielhaft die Existenz 
eines Nachbarvektors und präsentieren anschließend die Klassifikation und den 2- 
Nachbarschaftsgraphen des kompletten Geschlechtes des Barnes-Wall-Gitters. 

Eine Erzeugermatrix A für Aig entnehmen wir [CS92], S. 130. 
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1 / 



Wir betrachten den elften und zwölften Zeilenvektor aus der Matrix A und setzen 
V = 1, 1, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0). V ist ein Minimalvektor, also primitiv 

und insbesondere nicht in 2Ai6 enthalten. Das Skalarprodukt von v mit dem elften 

Basisvektor w = -^(2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0) ist gleich 1. Also ist w ^ 

(Ai6)t,, und damit ist gezeigt, daß (Ai6)i, / Aig. Da {v + 2w,v + 2w) = 4 + 4 + 16 = 
24 = (med 8), erhalten wir nach Lemma 2.5 mit Ai6(f + 2w) einen geraden 
Nachbarn von Aig. 
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2.7.2.1 Der Nachbarschaftsgraph 

Wir geben in Tabelle 2.6 den Nachbarschaftgraphen des Geschlechtes des Bames- 
Wall-Gitters an. Die Bezeichnung der Gitter folgt der Klassifikation in [SV94]. Die 
Vielfachheiten in der Inzidenzmatrix geben die Anzahl der Orbits an, die zu einer 
Nachbarbildung führen. 

2.7.3 Neuere Ergebnisse anderer Autoren 

2.7.3.1 ll-modulare Gitter in Dimension 12 

Ein extremales 11-modulares Gitter mitMinimum 8 in Dimension 12 wäre geringfügig 
dichter als das dichteste bekannte Gitter in dieser Dimension, das Coxeter-Todd- 
Gitter. Das Maß dieses Geschlechtes hegt bei ungefähr 15096.99 ([SH98], S. 748), 
aber die meisten Gitter besitzen relativ große Automorphismengruppen, so daß man 
aufgrund der Maßformel (siehe Gleichung 2.1 auf Seite 26)) von vielen hundert- 
tausend oder Millionen Gittern in diesem Geschlecht ausgehen kann. Eine vollständi- 
ge Klassifikation mit tn erscheint weder sinnvoll noch aussichtsreich. 

In [NV96] zeigen G. Nebe und B.B. Venkov, daß kein solches extremales Gitter ex- 
istiert. In dem Beweis, der sich im wesenthchen auf analytische Methoden und Siegel- 
Modulformen vom Grad 2 stützt, wird tn als Hilfsmittel benutzt, um ein Erzeugen- 
densystem des Vektorraums der vorkommenden Thetareihen zu konstruieren. 

2.7.3.2 Die Nachbarmethode im hermiteschen Faii 

In [Sch98] stellt A. Schiemann das Programm hn vor, daß die Nachbarmethode im 
hermiteschen Fall anwendet. Die erfolgreichen Klassifikationen sind in [Sch98] beschrieben. 

2.7.3.3 Extremale Gitter 

In dem Übersichtsartikel über extremale Gitter von R. Scharlau und R. Schulze-Pillot 
[SSP99] wurden die Ergebnisse mit tn und dem oben erwähnten hn berechnet. 

2.7.3.4 Diverses 

In der Doktorarbeit [Mis96] werden von M. Mischler einige Geschlechter mit tn 
bestimmt. 
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2.8 Verifikation der Ergebnisse 



Bei allen Klassifikationsproblemen stellt sich die Frage nach der Korrektheit und der 
Vollständigkeit des Ergebnisses. 

Das sogenannte Maß eines Geschlechts m{Q) kann einerseits mit analytischen Meth- 
oden hergeleitet werden (s. [CS88]), genügt andererseits nach Siegels Maßformel der 
Identität 



wobei die Summe über die Menge aller mit tn in ^ gefundenen Isomorphieklassen 
[L] gebildet wird. 

Diese Testfunktion erlaubt analog einer Quersumme in den oben zitierten Beispie- 
len eine sehr verläßliche Überprüfung der berechneten Klassifikationen. Die Maß- 
formel ist sensitiv gegenüber fehlenden Isomorphieklassen, einer fehlerhaften Berech- 
nung der (Ordnung der) Automorphismengruppen und inkorrekten Isometrietests. Eine 
fehlerhafte Berechnung eines Geschlechtes, die von der Maßformel nicht aufgedeckt 
wird, ist nur möglich bei gleichzeitigem Auftreten von zwei Fehlem, die sich bei der 
Berechnung der Maßformel gegenseitig aufheben. Die an diesem Test beteiligten Pro- 
grammteile sind die Generierung der Nachbarn, die Berechnung ihrer Automorphis- 
mengruppen und Isometrietests gegen die bisher gefundenen Gitter. Arbeiten zwei 
dieser drei Teile korrekt, so wird ein Fehler im dritten durch die Maßformel erkannt. 
Die Algorithmik des Isometrietests und der Berechnung der Automorphismengruppe 
sind sich in großen Programmteilen ähnlich, weil der Isometrietest mathematisch ver- 
gleichbar mit dem Auffinden eines nichttrivialen Automorphismus ist. In diesen Pro- 
grammteilen befinden sich aber eine ganze Anzahl von Plausibihtätstests, die Inkon- 
sistenzen entdecken können. 

Die wenigen bereits bekannten Klassifikationen (s.o.) konnten sowohl für unger- 
ade als auch für gerade Determinanten erfolgreich reproduziert werden. Die promi- 
nente Klassifikation der Niemeier-Gitter [Nie73] konnte mit tn nicht durchgeführt 
werden. Der Grund liegt in der algorithmischen Komplexität der Reduktion von Git- 
terbasen und in dem Umstand, daß die Automorphismengruppen in autom als Per- 
mutationsgruppen auf der Menge der kurzen Vektoren realisiert sind. Die Basis eines 
Nachbargitters, das mit dem oben beschriebenen Algorithmus konstruiert wird, en- 
thält zumeist einige längere Vektoren. In der Routine autom müssen alle Vektoren 
des Gitters bis zu der Länge des längsten Basisvektors berechnet werden. In der Di- 
mension 24 ist eine vollständige Aufzählung der kurzen Vektoren praktisch nur für 
Vektoren der Länge 2 und 4 möglich. Der LLL-Algorithmus reduziert aber nicht alle 
erzeugten Gitterbasen auf eine Basis, die nur Vektoren der Länge kleiner oder gleich 
4 enthalten. So bleibt die Klassifikation von Niemeier ein Meilenstein in der Theorie 




(2.1) 
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der ganzzahligen quadratischen Formen, für die es auch heute kein algorithmisches 
Äquivalent gibt. 
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Input: Eine Basis {61, • • • , bn} von L und v = {vi, - ■ ■ , u„) E L\2L mit (t;, f ) G 
4Z. 

Output: Eine Basis {6'^, • • • , 6^} des geraden Nachbarn von L bezüglich v, falls 

dieser existiert, ansonsten eine Fehlermeldung 

if {v, v) = 2 (mod 8) oder {v, v) = 6 (med 8) then 

retum FAILURE 
end if 

// 1. Initialisierung. 

if (v,v) = i (mod 8) then 

wähle i mit {v + bi,v + bi) = (mod 8) 

V *^ V + bi 
end if 

//2. Normalisiere v. 

k <— min{z | Vi ^ 2Z} 

wähle X mit x • = 1 (mod 8) 

V ^ X ■ V 

//3. Konstruktion der Basisvektoren. 
wähle m ^ k mit (i;, bm) ^ 2Z 
6^ < 2^) 26j„ 
for i = 1, • • • , n, fc 7^ i / m do 

if (6i,t;) G 2Zthen 
h'i^bi 

eise 

fe^ ^ 6j + 6m 

end if{// also {h\,v) £ 2Z } 
end for 
retum 6' 

Algorithmus 1: Eine Basis des geraden Nachbarn von L bzgl. v 
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3.1 Überblick 

Im Mittelpunkt dieses Kapitels steht ein Algorithmus nach einer geometrischen Idee 
M. Knesers zur Zerlegung eines Gitters in eine orthogonale Summe von unzerlegbaren 
Untergittem. Eine vereinfachte Variante dieser Methode ermögücht eine effiziente 
Konstruktion einer Gitterbasis aus einem großen Erzeugendensystem. Eine frühere 
Version dieses Kapitels ist in einen gemeinsamen Artikel [HV98] mit Frank Vallentin 
eingeflossen. 

In diesem Kapitel ist es zweckmäßig, ein Gitter L als Untermodul im n-dimensiona- 
len euklidischen Raum E versehen mit dem Standardskalarprodukt (., .) zu betrachten 
und den Aspekt der quadratischen Form als beschreibende Invariante eines Gitters in 
den Hintergrund treten zu lassen. Das Gitter L soll immer von vollem Rang sein, d.h., 
dim(L)]K = n. Im folgenden bezeichnen wir hier mit der Norm eines Vektors x e E 
seine euklidische Norm ||x|| = ^/{x,x). 

3.1.1 Das Zerlegungsproblem 

Bei vielen algorithmischen Konstruktionen von Gittern hat man nicht genügend Kon- 
trolle über die Eigenschaften und Invarianten des zu konstruierenden Gitters. Ins- 
besondere ist es für Klassifikationszwecke interessant, ob ein Gitter die orthogonale 
Summe zweier niederdimensionalerer Gitter ist. In dieser Situation benötigen wir 
einen effizienten Algorithmus, der ein Gitter in eine orthogonale Summe von Teilgit- 
tern zerlegt. 

Definition 3.1. Ein nichttriviales Gitter L heißt zerlegbar, falls es (echte) Untergitter 
Ü, L" C L gibt mit L = L' + L" und (Li, L2) = 0, andernfalls unzerlegbar. 

Bemerkung 3.2. Wir schreiben L = L' ® L" , falls L in L' und L" zerlegbar ist. 
In diesem Kapitel werden wir das Symbol © unter Strapazierung der Notation aus- 
schliesslich für eine innere direkte, orthogonale Summe verwenden. 

Die Existenz einer Zerlegung in unzerlegbare Untergitter folgt unmittelbar durch 
Induktion über den Rang von L. In einem nicht konstruktiven Beweis zeigte Eichler 
1952 in [Eic52] die Existenz einer bis auf die Reihenfolge der Summanden eindeuti- 
gen Zerlegung in unzerlegbare Untergitter. M. Kneser zeigte kurz darauf in einem 
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elegantem Beweis [Kne54], daß diese Zerlegung effektiv berechenbar ist. Mittels 
eines Erathostenessiebes werden Gittervektoren eliminiert, die nicht in einem unzer- 
legbaren Untergitter liegen. Damit erhält man ein Erzeugendensystem für die unzer- 
legbaren Untergitter. In Abschnitt 3.3 geben wir eine Version von Knesers Beweis für 
Erzeugendensysteme an, welche die benutzten geometrischen Methoden illustriert. 

Bei der Umsetzung von Knesers Beweis in einen Algorithmus stößt man auf zwei 
Probleme. Das Aussieben der zerlegbaren Vektoren ist ein aufwendiger Vorgang, 
dessen Laufzeit im wesentlichen quadratisch von der Größe des Erzeugendensystems 
abhängt. Des weiteren bildet die unzerlegbaren Vektoren im allgemeinen selbst ein 
sehr großes Erzeugendensystem für die unzerlegbaren Untergitter. Die Berechnung 
einer Gitterbasis aus einem Erzeugendensystem ist zwar ein algorithmisch gut ver- 
standenes Problem, in der Praxis jedoch für große Erzeugendensysteme aufwendig. 
Die bekannten Algorithmen wie der MLLL und die Bestimmung der Hermite-Normalform 
sind für Erzeugendensysteme konstruiert, die nicht wesentlich größer als der Rang des 
erzeugten Gitters sind. In der oben beschriebenen Situation führt dies in der Regel zu 
einem relativ schlechten Laufzeitverhalten. 

In Abschnitt 3.4 wird ein schnellerer Algorithmus für das Zerlegungsproblem als 
die Implementation von Knesers Methode angegeben; insbesondere erfolgt das Aussieben 
der erzeugenden Vektoren in linearer Laufzeit. 

3.1.2 Minimale Erzeugendensysteme 

Das zentrale Ergebnis in diesem Kapitel ist ein Überdeckungssatz für Erzeugenden- 
systeme, der eine obere Schranke für die Größe eines minimalen Erzeugendensystems 
aus Vektoren mit einer beschränkten Länge angibt. 

In Gittern mit nicht zu kleinem Minimum M (betrachtet im Verhältnis zur Determi- 
nante) besteht jedes minimale Erzeugendensystem nur aus wenig mehr als n Vektoren. 
Der Überdeckungssatz 3.14 quantifiziert dieses Ergebnis genauer: Ist 5 C L ein min- 
imales Erzeugendensystem von L mit Vektoren, die nicht länger als ß G M sind, dann 
ist|5Kn + log2(n!(f r). 

Die algorithmische Hürde dieser Aufgabe rührt vom Fehlen eines Basisergänzungssatzes 
und eines Steinitz'schen Basisaustauschsatzes für die Modulsituation her. Man kann 
sich dies leicht an einem Erzeugendensystem veranschaulichen, in dem kein Vektor 
im Erzeugnis der anderen hegt, siehe etwa das Beispiel 3.21 auf Seite 44. 

3.1.3 Der l\/ILLL bei großen Erzeugendensystemen 

Die in Kapitel 3.4 entwickelten Ideen für eine effiziente Lösung des Zerlegungsprob- 
lems werden in Kapitel 3.5 eingesetzt, um einen effizienten Algorithmus zu kon- 
struieren, der aus einem großen Erzeugendensystem eine Gitterbasis berechnet. Die 
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Methode besteht aus einem Meta-Algorithmus, der als Kern bekannte Routinen wie 
den MLLL oder die Berechnung der Hermite-Normalform benutzt. Diese Kemalgo- 
rithmen werden wiederholt aufgerufen, jedoch ausschüeßlich für minimale, also sehr 
kleine Erzeugendensysteme. Für den Meta-Algorithmus können bei großen Erzeugen- 
densystemen schärfere asymptotische Schranken gezeigt werden als für die zugrun- 
deliegenden Kernalgorithmen. Die praktische Relevanz erkannt man an den Laufzeit- 
gewinnen, die in den graphischen Auswertungen zweier Experimente (s. S. 45, Abb. 3.3 
und Abb. 3.4) veranschaulicht werden. 

3.2 Unzerlegbare Vektoren 

Bevor wir uns der konkreten Berechnung einer Zerlegung von L in unzerlegbare Un- 
tergitter zuwenden, betrachten wir quaütative Kriterien für die Zerlegbarkeit bzw. Un- 
zerlegbarkeit von Vektoren. Für ein Gitter mit der Zerlegung L = Li . . . © sollen 
unzerlegbare Vektoren zwei Eigenschaften genügen: 

i) Wenn x e Li, y e Lj,i ^ j unzerlegbare Vektoren sind, dann sei x -|- y G 
Li © Lj zerlegbar. 

ii) Die Menge der unzerlegbaren Vektoren in Lj, 1 ^ z ^ r erzeugt Lj. 

Der Begriff der Unzerlegbarkeit von Vektoren wurde zuerst von M. Kneser definiert 
[Kne54] und später von O'Meara [O'M80] verfeinert.^ 

Definition 3.3. Sei v e L. 

i) V heißt orthogonal zerlegbar, /«//.v Vektoren x,y £ L\0 existieren mitv = x + y 
und {x, y) = 0, andernfalls orthogonal unzerlegbar. 

ii) V heißt spitz zerlegbar, /a//* Vektoren x,y & L\0 existieren mit v = x + y und 
(x, y) ^ 0, andernfalls spitz unzerlegbar. 

iii) V heißt linear zerlegbar oder einfach zerlegbar, falls Vektoren x,y £ L existieren 
mit V = X + y und \ \v\ \ > \ \x\ \ ^ \ \y\\, andernfalls (linear) unzerlegbar. 

Lemma 3.4 ([O'M80]). Sei v e L, dann gih: 

V unzerlegbar =^ v spitz unzerlegbar =^ v orthogonal unzerlegbar 

Beweis. Sei v orthogonal zerlegbar, dann ist v offenbar auch spitz zerlegbar. Sei v 
nun spitz zerlegbar, d.h., es gibt x,y e L\0 mit v = x + y und (x, y) ^ 0. Dann 
folgt aus 1 1^1 1 = ^/x^'^\^'2(^c^yy^^y^ > \\x\\, daß v zerlegbar ist. □ 

'wir danken M. Kneser für seinen freundlichen Hinweis auf die Arbeit von O'Meara. 
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O'Meara gibt eine Charakterisierung orthogonal und spitz zerlegbarer Vektoren an. 
Satz 3.5 ([O'M80] 3.6, 3.7). Sei v e L. 

i) V ist genau dann orthogonal zerlegbar, wenn ein w & L existiert mit w ^ 

±v, V — w G 2L und \\v\\ = \\w\\. 

ii) V ist genau dann spitz zerlegbar, wenn ein w E L existiert mit w ^ ±v, v — wE 

2L und \ \v\\ ^ \\w\\. 

Beweis. Wir folgen O'Mearas Beweis und zeigen zunächst Aussage ii). Sei v e L. 
V spitz zerlegbar 

(Nun sei w := v — 2y bzw. y ■= ^{v — w)). 

3w e L \ {±v} {^{v - w), ^{v + w)) ^ Ound v - w e 2L 

& L \ {iv} {v, v) — {w, w) und v — w e 2L 
3w & L \ I |i;| I ^ I \w\ I und v — w e 2L 



Aussage i) erhalten wir, wenn wir im vorangehenden Beweis „spitz" durch „orthog- 
onal" und alle „^" durch „=" ersetzen. □ 




Für v,y G L mit {y,v — y) > ex. ein w e L mit | |w| | < \ \v\ \ und v — w G 2L. 
Abbildung 3.1: Illustration zum Beweis von Satz 3.5 ii) „^" 
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Damit hat O'Meara eine bemerkenswerte Eigenschaft von voronoirelevanten Vek- 
toren bewiesen, allerdings ohne den Zusammenhang zur Voronoitheorie zu erwäh- 
nen.-^ Unseres Wissens ist dieser simple Zusammenhang bis heute nicht formuliert 
worden. 

KoroUar 3.6. Ein Vektor ^ x E L ist genau dann spitz unzerlegbar, wenn er 
voronoirelevant ist. 

Beweis. Nach Satz 3.5 ist x genau dann spitz unzerlegbar, wenn ±x die einzigen 
kürzesten Vektoren in der Nebenklasse x + 2L sind. Die Behauptung folgt nun un- 
mittelbar aus einem Kriterium von Voronoi, hier Satz 4. 12, S. 52. □ 

Aus Satz 3.5 folgt ebenfalls, daß es höchstens 2 • (2"^ — 1) spitz unzerlegbare 
(und damit hnear zerlegbare) Vektoren gibt ([O'M80] 3.17, 3.18). Eine einfache ge- 
ometrische Überlegung zeigt bereits, daß lange Vektoren immer zerlegbar sind. 

Lemma 3.7. Sei R der Überdeckungsradius^ von L und v & L mit \ \v\\ > 2R. Dann 
ist V zerlegbar 

Beweis I:. Seiw e L mit — ^ R. Dann gilt 

1 ,, ,,1 ,, „ 1 

||u;|| ^ Ww — -v\ \ + 112^11 ^ -f^ + 2 ll"^!! ll^ll 



und 



1 ,, ,, 1 ,, „ 1 

Kf^ — ^ \\w — + v\\ ^ R-\ — b < b 



Also ist V zerlegbar, weil v die Summe der beiden echt kürzeren Vektoren w und v — w 
ist. □ 

Mit einem Argument aus der Voronoitheorie folgt die Behauptung unmittelbar. 

Beweis II:. (F. Vallentin) Ist > 2R, so folgt aus ^ Vb(-L), daß v nicht 
voronoirelevant ist und mit Satz 3.5 zerlegbar sein muß. □ 

Bemerkung 3.8. Der Begriff der orthogonalen Zerlegbarkeit von Vektoren wurde von 
M. Kneser analog zum Begriff der Zerlegbarkeit von Gittern definiert [Kne54]. In 
generischen Gittern sind alle Vektoren orthogonal unzerlegbar. Ein einfaches Beispiel 
hierfür ist das von den Vektoren (1, 0), {9, 1) mit ^ Q erzeugte Gitter. Insbeson- 
dere sind in einem eindimensionalen Gitter alle nichttrivialen Vektoren orthogonal 
unzerlegbar. O'Mearas Definition der linearen Zerlegbarkeit von Vektoren erscheint 
uns intuitiver. 



^Siehe die Definition 4. 10 auf Seite 52. 
Der Überdeckungsradius ist der m; 
in L, siehe Definition 4.3, S. 49. 



^Der Überdeckungsradius ist der maximal mögliche Abstand eines Punktes in E von einem Gitterpunkt 
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V 




Abbildung 3.2: Lemma 3.7: Lange Vektoren sind zerlegbar 

3.3 Knesers Methode 

Sei B £ R, so daß S := {^^i, • • • ,Vs} = € L \ | ||f || ^ B} ein Erzeugen- 
densystem von L ist. Auf der Menge der unzerlegbaren Vektoren I in S definieren 
wir den Orthogonaütätsgraph T = (/, E) mit E = {{v,w} e I x I \ {v,w) ^ 0}. 
Wir zerlegen T in seine Zusammenhangskomponenten Fj = (7j, 1 ^ i ^ r und 
setzen = (/j)z. 

Satz 3.9 ([Kne54]). Die Gitter Li,l ^ i ^ r sind unzerlegbar und die orthogonale 
Zerlegung 

L = Li®...®Lr (3.1) 
ist eindeutig bis auf die Reihenfolge der Summanden. 

Beweis. Seiv £ S \ und \ \v\ \ minimal mit dieser Eigenschaft. Dann ist v wegen 
V ^ I zerlegbar und es gibt x,y e S von echt kürzerer Norm als v mit v = x + y. 
Mindestens einer dieser beiden Vektoren ist nicht in {I)z, sonst wäre auch v G 
Dies widerspricht jedoch der Minimalität von II II, somit ist 5 C gezeigt. Daraus 
folgt die erste Behauptung L = Yli=i ^i- 

(Ij, Ij) = für i 7^ j impliziert {Li, Lj) = 0. Daraus folgt mit der Definitheit von 
(., .) und der vorhergehenden Aussage, daß 

L = Li® Lr. 

Seien L' , L" C Li für ein beliebiges, aber festes i mit Li = L' ® L" . Für v £ Ii 
gilt entweder v £ L' oder v £ L", ansonsten wäre v zerlegbar. Ohne Beschränkung 
der Allgemeinheit sei v G L'. Wegen {L\L") = sind L' DT und L" D F zwei 
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nicht zusammenhängende Teilmengen von F. Also ist L" = und L' = L^. Damit ist 
gezeigt, daß Li unzerlegbar ist. 

Da jedes Lj bereits durch einen in ihm enthaltenen unzerlegbaren Vektor eindeutig 
bestimmt ist, folgt die Eindeutigkeit der Zerlegung bis auf Reihenfolge der Summan- 
den. □ 

Wir beschreiben nun Knesers Methode zur Konstruktion von T : 

Initialisierung Sei {61, • • • , 6„} eine Basis von L. Setze B = maxj und berechne 

S = {veL\0\ \\v\\ ^B}. 

Konstruiere / Für jedes Paar u,v & S mit u + v & S markiere u + v als zerlegbar. 
Sei I die Menge aller nicht als zerlegbar markierten Vektoren. 

Konstruiere E Füge {u, v} für alle u,v e I inE ein, falls (u, v) / 0. 

Bei diesem Verfahren müssen wir jedes Paar {u, v} für alle u,v E S untersuchen. 
Diese (2) G O(s^) Operationen sind für große s sehr aufwendig. 



3.4 Der inkrementelle Algorithmus 

Um für ein Lj in Gleichung 3. 1 eine Basis effizient zu bestimmen, ist es nicht notwendig, 
alle unzerlegbaren Vektoren in Li zu berechnen. Es genügt, ein Erzeugendensystem 
in den Lj aufzufinden, welches einerseits aus relativ wenigen Vektoren besteht und 
andererseits erlaubt, schnell zu verifizieren, daß sich alle unzerlegbaren Vektoren in 
den daraus erzeugten Gittern befinden. Dies motiviert den folgenden Algorithmus 2. 
Die Korrektheit dieses Algorithmus zeigt 

Satz 3.10. Sei L ein n-dimensionales Gitter auf E, B E so daß die Menge S = 
{vi, ■ ■ ■ ,Vs} := {v G L \ \\v\\ ^ B} ein Erzeugendensystem von L enthält. Dann 
berechnet Algorithmus 2 die Zerlegung von L in unzerlegbare Untergitter Li, d. h., 
L = Li® . . .® Lr. 

Für ein gegebenes System von paarweise orthogonalen Untergittem Li © • • • © 
sei TTj die orthogonale Projektion auf den R- Vektorraum, der von Lj aufgespannt wird. 
/ und r seien wie oben definiert. 

Beweis. Sei Li ® . . . ® Lk mit den folgenden Eigenschaften bereits berechnet: 

i) L = (.s)z + Eti^i, 

ü) Li ist unzerlegbar für 1 ^ z ^ A;, 
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1: Input: ß G M mit einem Erzeugendensystem S* = {f G L \ | 1 1^1 1 ^ B} von L. 

2: Output: Unzerlegbare Untergitter Lj mit L = Li © . . . L^. 

3: // 1. Initialisierung. 

4: Wähle v G S mit | \v\ \ minimal. 

5: S\{v}, k^l,Lk^ 

6: // 2. Füge sukzessiv Vektoren aus S in die Li ein. 

7: whUe S* / do 

8: Wähle v E S mit minimalem 1 1^1 1. 

9: S ^ S\{v}. 
10: if -ü Li © . . . © Lfe then 
11: // V is „ergänzend", insbesondere unzerlegbar. 
12: J ^ {i e {!,■■■ ,k}\Tri{v)^0} 
13: L^Zz; + 0.gjLi 
14: // Ordne die Liste der Li neu. 

15: {Li,--- ,Lfe_|j|} ^ {Lj I j ^ J}, Lfe_|j|+i^L, k^k-\J\ + l 
16: end if 

17: end while 

Algorithmus 2: Zerlegung des Gitters L 

iii) (Li, Lj) = für i / j. 

Diese Eigenschaften sind nach der Initialisierung in Zeile 5 von Algorithmus 2 erfüllt. 

Wir zeigen nun die Schleifeninvarianz (Zeile 7-16) dieser Eigenschaften. Sei v ^ S 
ein Vektor mit minimaler Norm in S wie in Zeile 8 gewählt. Ist G Li © . . . © L^^.) 
so können wir v verwerfen und die Eigenschaften i) - iii) bleiben erhalten. Sei nun 
u ^ Li © . . . © Ljfc. Einen solchen Vektor neimen wir „ergänzend"^. 

Zu i): V wird hinzugefügt und es gilt 

k 

L={S\{v}U + {I.v + Y,Li). 

i=l 

Insbesondere gilt L = Li © . . . © Lj. nachdem alle Vektoren in S abgearbeitet 
sind. 

Zu ii): V ist unzerlegbar, sonst gäbe es x,y G S von echt kürzerer Norm als v mit 
v = X + y. Nach Voraussetzung wären x,?/ jedoch schon in einer vorherigen 

'*Die Begriffsbildung ist lediglich lokal, d.h. nur im jeweiligen Schleifendurchlauf sinnvoll, weil diese 
Eigenschaft von der Menge der bereits abgearbeiteten Vektoren in 5 abhängt und nicht eine Invari- 
ante des Vektors v ist. 
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Schleife abgearbeitet worden, woraus x + y = vGLi(B-..(BLk folgen würde. 
(Insbesondere folgt aus der orthogonalen Summe v = {v — 'Ki{v)) + ni{v), 1 ^ 
i ^ k, daß es keine nichttriviale Projektion von v auf {Li)Q gibt, welche in Li 
hegt.) Wir setzen J = {j E {!, - ■ ■ ,k} \ T^j{v) ^ 0} und wählen Vektoren 
Vj G Lj,j € J mit {vj,v) ^ 0. Dann ist der Teilgraph {{vj. v)}j(zj C V ein 
Stern und damit zusammenhängend, also ist 0jg j Lj + "Lv unzerlegbar. 

Zu iii): 0jg/y Li © (0jgj Lj + 'Lv) ist eine orthogonale Zerlegung, weil 'Kiiy) = 
für i ^ J. 

Damit ist die Korrektheit von Algorithmus 2 gezeigt. □ 
3.4.1 Datenstrukturen und Laufzeitanalyse 

In diesem Abschnitt beschreiben wir für die Reahsierung von Algorithmus 2 geeignete 
Datenstrukturen. Wir werden in den Sätzen 3.17 und 3.19 Laufzeitanalysen für die 
hier entwickelten Algorithmen aufgrund dieser Datenstrukturen vornehmen. Im fol- 
genden fassen wir E als Vektorraum aus Spaltenvektoren auf. 

Eine zentrale Rolle in diesem Algorithmus spielen die orthogonalen Zerlegungen 
-Li © . . . -Ljfc, die in jedem Schleifendurchlauf konstruiert werden. Jedes Teilgitter 
Lj wird durch eine Gitterbasis vn, - ■ ■ , Ui^j C E repräsentiert. Des weiteren benöti- 
gen wir eine (möglicherweise leere) Vektorraumbasis wi, - ■ ■ , Wn-i des orthogonalen 
Komplements (Li, • • • , L/j)^, l = n — Yli=i ''i- Nun setzen wir 

A = (yu ■ ■ ■ Vkvk Wi ■■■ Wn-l) ^ e ajtnxn W- 

Für jedes x G E gilt v G Li©...©Lfc genau dann, wenn die ersten n — l Ko- 
effizienten von Av in Z liegen und die resthchen l gleich sind. Wir können somit 
in Zeile 10 des Algorithmus durch eine Matrix- Vektor-Multiplikation testen, ob der 
Vektor x in dem bereits konstruierten Untergitter liegt. Ist das der Fall, so übersprin- 
gen wir ihn und gehen zum nächsten Schleifendurchlauf. Andernfalls berechnen wir 
wie folgt das neue Gitter L mit einer seiner Basen, ordnen die Liste unserer unzerleg- 
baren Teilgitter Lj neu und berechnen die Indexmenge / und dann die Matrix A in der 
folgenden Weise neu: A ist eine Diagonalmatrix mit den Blöcken [vn ■ ■ ■ Vir^) ^ 
auf der Hauptdiagonalen. Für jeden der ersten l Koeffizienten von Av, welcher nicht 
in Z liegt, suchen wir den entsprechenden Block auf (etwa j) und fügen j in J ein. 
Anschließend konstruieren wir eine Basis von L = Ylj^J einer Variante^ des 

bekannten LLL- Algorithmus [LLL82]. Die Berechnung einer Basis von L kann ef- 
fizient erfolgen, weil L ein Erzeugendensystem {vji | j € J, 1 ^ i ^ r,} U {v} aus 

^Beispielsweise mit einer der in [Poh87] oder [BP89] beschriebenen Modifikationen. Siehe Ab- 
schnitt 3.5 für eine ausführliche Diskussion der Alternativen. 
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höchstens ^j^jrj ^ n + 1 Vektoren besitzt. Nun ordnen wir die Liste der Gitter 
L,Li,i ^ J neu und indizieren sie von 1 bis — | J| + 1. SchUeßlich berechnen wir 
A mittels Gaußtransformationen neu. 

Die Effizienz dieses Algorithmus beruht auf der Tatsache, daß sich die meisten 
V ^ S in einer bereits konstruierten Zerlegung Li © ... © befinden und die 
aufwendige Berechnung des Gitters M und seiner Basis im Verhältnis zur Mächtigkeit 
von S nur selten notwendig ist. Im folgenden werden wir diese Behauptung genauer 
quantifizieren. Dazu benutzen wir hauptsächüch Minkowskis zweiten Fundamental- 
satz. 



3.4.1.1 Ketten von Untergittern 

Satz 3.11 ([Min96]). Sei L ein Gitter auf dem n-dimensionalen euklidischen Raum 
E mit den sukzessiven Minima Ai(L), • • • , Xn{L) und = {x G E \ \\x\\ ^ 1} der 
n-dimensionale Ball mit Radius 1. Dann gilt 

2" 

—VdetL ^ Ai(L) • . . . • XJL) ■ volBn ^ 2"VdetL. (3.2) 
n! 

Im Zentrum dieses Abschnittes steht der folgende 

Satz 3.12. Sei L ein Gitter mit Minimum M auf dem n-dimensionalen euklidischen 
Raum V. Das Untergitter L' C. L sei ebenso wie L von Vektoren erzeugt, die alle nicht 
länger als B E sind. Weiter gebe es eine endliche Kette von Untergittern Lj C L 
mit 

L' = Li C ■■■ C Lt = L (3.3) 
Dann gilt für die Länge t der Kette: 
i) Ist rankL' = rankL, dann ist 

(B \ 
m) 



//) Wenn alle Li aus Vektoren nicht länger als B erzeugt werden, dann ist 

t + log2(n! 



M J 



n 
I ) 



Beweis. Zu i): Sei Sj = [L : Li] der Index von Lj in L, 1 ^ z ^ t. Dann gilt nach der 
Determinanten-Index-Formel sf = [L : Lif = ^^j;, daß 

VdetL'/ det L = [L : L'] = si > • • • > st = [L : L] = 1. 
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Sj+i ist wegen ,Si+i • [Li+i : Li] = ,Sj ein Teiler von s«. Also ist die maximale Anzahl 
von Gittern in der Kette (3.3) durch die längstmöghche Teilerkette von det L' / det L 
nach oben beschränkt. Der kleinstmögüche Teiler ist jeweils 2, also erhalten wir für t 
eine obere Schranke von logg y^detL'/ detL. 
Nach Satz 3.11 gilt für L' 

VditL^ ^ - voli?„ • Ai(L') • . . . • A„(L') < — volB^B^ 

und für L 

VditL > ^voli?„ • Ai(L) • . . . • A„(L) ^ ^YolBnM^. 

Daraus folgt 

Zu ii) Wir nehmen ohne Beschränkung der Allgemeinheit an, daß rankLj+i — 
rankLj ^ 1 für alle 1 ^ i < n ist. Ansonsten verfeinern wir die Kette und fü- 
gen neue Gitter hinzu, bis sich zwei in der Kette benachbarte Gitter jeweils im Rang 
um höchstens 1 unterscheiden. Dann zeigen wir die Behauptung für die verfeinerte, 
längere Kette. 

Nun werden die Gitter umnumeriert, um sie mit ihrem Rang zu indizieren. Die 
Kette bestehe nun aus Gittern Lij C L, 1 ^ j ^ vom Rank i, erzeugt von 
Vektoren nicht länger als ß G M, so daß 

Loi C Lu C Li2 C • • • C Ln,mn-1 C Lnm„ = L. (3.4) 

Für 1 ^ z ^ n sei € Lu \ (ij-i,mi_i)R mit ||2;j|| ^ B. Ein solches Zi existiert, 
weil alle Gitter Lij aus Vektoren nicht länger als B erzeugt werden. Wir setzen F = 
(zi, • • • , Zn)z- Die Menge {zi, • • • , Zn} ist nach Konstruktion M— linear unabhängig, 
also ist rank F = n. Nun addieren wir F zu jedem Gitter in der Kette (3.4) hinzu. 
Damit erhalten wir die folgende Kette von Gittern, die alle den Rang n besitzen. 

{F =)Loi + F C . . . C Li,- + F C . . . C Lnmr. + F{= L). (3.5) 

Man verifiziert leicht, daß L^j + F ^ + j < fi für alle 1 ^ i ^ n gilt. Dann 

kann in der Kette (3.5) Gleichheit nur an den Positionen Li_i y.. ^ C Lu +F, 1 < 
i < n, also höchstens n-mal gelten. Wir wenden die in Teil i) gezeigte Aussage auf 
die Kette der verschiedenen Gitter in (3.5) an und haben damit gezeigt, daß die Kette 
(3.4) und damit auch die Kette in Behauptung ii) nicht mehr als n + log2(n!(;^)") 
Gitter enthält. □ 
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3.4.1.2 Ein Überdeckungssatz für Erzeugendensysteme 

Wir fassen diese Ergebnisse in einem „Überdeckungssatz" zusammen, der die Exis- 
tenz von kleinen Erzeugendensystemen garantiert. 

Definition 3.13. Ein Erzeugendensystem S eines Gitters heißt minimal, wenn für alle 
S' gilt 

S' cS^ {S')z c {S)z. 

Satz 3.14, Sei L ein Gitter mit Minimum M auf dem n-dimensionalen euklidischen 
Raum E. Sei S ein minimales Erzeugendensystem von L mit Vektoren, die nicht länger 
als B eM. sind. Dann ist 

\S\^n + log,{n\{^r). 

Beweis. Sei S = {si, • • • , sj}. Wegen der Minimaütät von S gilt (si, • • ■ , Si)z 7^ 
(si, • • • , für alle 1 ^ i < t. Die Behauptung folgt aus der Anwendung von 

Satz 3.12 ii) auf die Kette 

(Sl)z C (Sl, • • • , Si)z C (si, • • • , st)z = {S)z- 

□ 

3.4.1.3 Laufzeitanalyse von Aigoritlinrius 2 

Eine asymptotische Schranke für die Anzahl der „teuren" Updateoperationen in Algo- 
rithmus 2, bei denen die Basis von L und die Matrix A, wie auf Seite 37 beschrieben, 
neu berechnet werden müssen, erhalten wir in 

Korollar 3.15. Die Anzahl der Updateoperationen in Algorithmus 2 für die Zerlegung 

eines Gitters L ist in 0(nlog 

Beweis. In Algorithmus 2 werden sukzessive Zerlegungen von Teilgittern Li © ... © 
Lk ^ L konstruiert. Diese Teilgitter erfüllen die Voraussetzungen von Satz 3.12 ii). 
Damit ist die Anzahl der ergänzenden Vektoren, die eine Updateoperation notwendig 
machen, nicht größer als n + log2(n!(;^)"). Auf diese Schranke wenden wir die 
Stirhng'sche Formel n! \/27rn(^)" an und erhalten eine obere Abschätzung, die in 
0(n log ^) liegt. □ 

Ein wichtiger Schritt ist die Konstruktion des Gitters L im Algorithmus 2, Zeile 13. 
Wir benötigen hierfür ein Verfahren, das aus einem Erzeugendensystem von höch- 
stens n + 1 Vektoren eine Basis von L effizient berechnet. Weiterhin sollten diese 
Basisvektoren in irgendeinem Sinn reduziert, d.h. heißt möglichst kurz sein. Ohne 
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diese Reduktion kann es bei der Iteration von Basisberechnungen zu einer Koeffizien- 
tenexplosion kommen, welche die Laufzeit der späteren Basisberechnungen außeror- 
dentlich verschlechtert. Dieses Verhalten ist beispielsweise bei der Berechnung der 
Hermite-Normalform problematisch (siehe das Beispiel von Hafner und McCurley in 
[HM89]). Wir ziehen daher die Verwendung einer Variante des LLL-Algorithmus 
vor, der bereits während der Basisberechnung reduziert. Es sind exakte Schranken 
für die Länge von LLL-reduzierten Basisvektoren in Beziehung zu den sukzessiven 
Minima des Gitters bekannt. Für die maximale Laufzeit des LLL-Algorithmus ist bei 
ganzzahligen Gittern eine Abschätzung bekannt. 

Lemma 3.16. Sei L ein ganzzahliges Gitter auf dem n-dimensionalen euklidischen 

Raum E, welches von n + 1 Vektoren, alle kürzer als i? G M, erzeugt wird. Dann 
kann eine (LLL-reduzierte) Basis {5i, • • • , 6„} von L in 0{n^ log B) arithmetischen 
Operationen berechnet werden und es gilt für 1 ^ j ^ n, daß 

\\hj\\ ^ 2("-i)/2 . A,(L) ^ 2('^-i)/2 . ^^^^y (3_6) 

Beweis. Nach [BP89], Theorem 3.1. läßt sich eine LLL-reduzierte Basis in 0(2n + 
2Y log B C 0{n^ log B) arithmetischen Operationen berechnen. Die Länge der Ba- 
sisvektoren ist nach [LLL82], Proposition 1.12 wie behauptet beschränkt. □ 

Nach diesen Vorbereitungen können wir die Laufzeit von Algorithmus 2 nach oben 

abschätzen. 

Satz 3.17. Sei L ein ganzzahliges Gitter auf dem n-dimensionalen euklidischen Raum 
E,B gW, so daß S := {vi , • • • ,Vs} = {v & V\0 \ \\v\\ ^ B} das Gitter L erzeugt. 
Dann können wir eine orthogonale Zerlegung L = Li ... von L in unzer- 
legbare Teilgitter Li,l ^ i ^ r in höchstens 0{nPlog^{nB) + sn?) arithmetischen 
Operationen berechnen. 

Beweis. Wir analysieren die Laufzeit von Algorithmus 2. Dazu benutzen wir die in 
Abschnitt 3.4.1 beschriebenen Datenstrukturen. 

Mittels Radix-Sort ([Knu74]) sortieren wir alle Vektoren in S ihrer Norm nach und 
fügen sie in eine Liste ein. Die Norm berechnen wir in 0(n) Operationen, so daß der 
gesamte Sortiervorgang in 0{sn) Operationen erfolgen kann. 

Seien nun in Zeile 7 bei Eintritt in die Schleife vu, • • • , vtr^,wi, • • • , Wn^i und 
A bereits per Induktion berechnet. In Zeile 8 wird ein v G S mit minimaler Norm 
entnommen und geprüft, ob v G Li ... L^, d.h., ob die ersten n — l Koef- 
fizienten von Av in Z liegen und die restlichen l Koeffizienten gleich sind. Dies 
kostet für alle v £ S zusammengenommen O(sn^) Operationen. Falls diese Bedin- 
gung erfüllt ist, dann liegt v in dem bereits erzeugten Untergitter und es wird kein 
Update der Datenstrukturen benötigt. Nach Satz 3.12 ii) werden für ganz S in keinem 
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Fall mehr als ^ log(n!(-p)"^) + n Updatevorgänge benötigt. Bei jedem Update wird 
eine LLL-reduzierte Gitterbasis aus einem Erzeugendensystem mit höchstens n + 1 
Vektoren berechnet. Diese Vektoren sind entweder aus S oder Output eines früheren 
LLLs, also LLL-reduziert. Damit ist ihre Norm durch 2("-i)/2A„(L) ^ 2("-i)/2ß 
nach oben beschränkt. Nach Lemma 3.16 können wir eine Basis des erzeugten Git- 
ters in höchstens 0{n* log(2("~^)/^ß)) = 0(n^ + log B) Operationen berechnen. 
Die Matrix A können wir beispielsweise mit dem Gaußverfahren (siehe [GVL96]) in 
O(n^) Operationen updaten. 

In der Summe erhalten wir damit für die Gesamtlauf zeit eine obere Schranke von 
0(n^ log{nB) + log{nB) log B + .sn^) arithmetischen Operationen. Daraus folgt 
die leicht schlechtere, aber in der Schreibweise übersichtlichere Schranke von 0(71^ log^(nß)+ 
sn^) arithmetischen Operationen. □ 

Bemerkung 3.18. Conway und Sloane konstruierten 1995 ein Gitter, das aus seinen 
Minimalvektoren erzeugt wird, aber keine Basis aus Minimalvektoren besitzt [CS95]. 
Damit ist gezeigt, daß es nicht möglich ist, das Anwachsen der Normen der Basisvek- 
toren in der neuen Basis gegenüber dem ursprünglichen Erzeugendensystem zu ver- 
hindern. Dennoch ist die Abschätzung über die Länge von LLL-reduzierten Vektoren 
in Gleichung (3.6) äußerst pessimistisch. In der Praxis findet der LLL-Algorithmus 
deutlich kürzere Basisvektoren, insbesondere wenn das Erzeugendensystem bereits 
aus relativ kurzen Vektoren besteht. 

3.5 Gitterbasen aus großen Erzeugendensystemen 

Ein in der Praxis wichtigeres und viel häufigeres Problem als die Zerlegung von Git- 
tern ist die Berechnung einer Gitterbasis aus einem großen Erzeugendensystem. Hi- 
erzu stellen wir die bekanntesten Algorithmen samt einiger Modifikationen zusam- 
men. Im folgenden sei L ein ganzzahliges Gitter mit Minimum M im n-dimensionalen 
euklidischen Raum E, das von s Vektoren erzeugt wird, die alle nicht länger als 
B eM. sind. Weiter nehmen wir an, daß s > n gilt: 

3.5.1 DerBP-LLL 

Die schon in Lemma 3.16 angesprochene Modifikation des LLL-Algorithmus von 
Buchmann und Pohst ([BP89]), im folgenden BP-LLL genannt, wird in der Praxis 
selten eingesetzt. Für theoretische Überlegungen ist dieses Verfahren wertvoll, weil 
dafür eine obere Abschätzung der Laufzeit von 0{{s + n)^logB) C 0{s^logB) 
arithmetischen Operationen existiert. 
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3.5.2 Der MLLL 

Die von Pohst in [Poh87] beschriebene Modifikation des LLL ist unseres Wissens 
der meistverwendete Algorithmus zur Basisbestimmung. Leider ist die theoretische 
Laufzeit des MLLL gegenwärtig nicht bekannt. In der Praxis ist er jedoch in seinem 
Laufzeitverhalten für kleine Erzeugendensysteme (d. h. s n) dem LLL ähnüch. 

3.5.3 Der serielle BP-LLL 

Anstatt den BP-LLL auf das gesamte Erzeugendensystem mit s Vektoren anzuwen- 
den, benutzen wir ihn sukzessiv für Erzeugendensysteme bestehend aus n + l Vek- 
toren. Die Länge der Vektoren in den dabei konstruierten Basen lassen sich mit 
Lemma 3.16 durch 2("-i)/2A„(L) ^ 2("-i)/2ß nach oben abschätzen. Damit ist 
die Laufzeit des seriellen BP-LLL in 0{s{2n + 2)^log(2("-i)/2ß)) C 0(sn''(n + 
logß)). 

3.5.4 Der serielle BP-LLL mit Aussieben 

Wir können den Algorithmus 2 vereinfachen, um lediglich eine Gitterbasis statt einer 
Zerlegung des Gitters zu berechnen. Der daraus entststehende Algorithmus 3 (s. S. 43) 
benutzt eine Unterfunktion construct_basis, die eine Basis aus höchstens n + 
1 Erzeugern berechnet. Um die Laufzeit mit Satz 3.19 abzuschätzen, benutzen wir 
hierfür den BP-LLL; man beachte aber auch hier die Bemerkung in 3.5.2. 



1 


Input: Ein Erzeugendensystem S = {vi, ■ ■ ■ ,Vs} von L mit \ \vi\\ ^ B. 


2 


Output: Eine Basis • • • , bn} von L. 


3 


1. //Initialisierung. 


4 


Wähle veS,S^ S\{v}, bi ^ v 


5 


2. //Füge Vektoren v aus S sukzessive hinzu. 


6 


while 5 / do 


7 


Wähle G 5 und setze S ^ S\{v}. 


8 


ifv^ {bi, ■ ■ ■ ,bk)z then 


9 


{6i,--- jfejfc}-«— construct_basis(6i, • • • ,bk,v) 


10 


end if 


11 


end while 


Algorithmus 3: Berechnung der Basis von L aus einem Erzeugendensystem S 



Satz 3.19. Sei L ein ganzzahliges Gitter auf dem n-dimensionalen euklidischen Raum 
E, das von s Vektoren, alle nicht länger als B & R, erzeugt wird. Dann können wir 
eine Basis von L in höchstens 0(n^ log^(nZ?) + sn^) arithmetischen Operationen 
berechnen. 
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Beweis. Wir benutzen für die Formulierung von Algorithmus 3 dieselben Datenstruk- 
turen wie sie in Abschnitt 3.4.1 beschrieben sind und verwenden den BP-LLL für die 
Unterfunktion construct_basis. Dann berechnet der folgende Algorithmus 3 
eine (LLL-reduzierte) Basis von L. Die Abschätzung der Laufzeit erfolgt vollständig 
analog zum Beweis von Satz 3.17. □ 

Bemerkung 3.20. Im Vergleich zum seriellen BP-LLL werden in Algorithmus 3 für 

große s die meisten LLLs durch eine einzige Matrix- Vektor-Multiplikation ersetzt. 
Dafür werden die im Vergleich zu s relativ wenige Updates mit der Berechnung der 
Inversen einer n x n-Matrix verteuert. 



3.6 Beispiele und Beobachtungen 

Beispiel 3.21. Sei E = M.^ und {pi, • • • ,pk} eine Menge von paarweise verschiede- 
nen Primzahlen, deren Produkt gleich t ist. Dann folgt mit dem chinesischen Reste- 
satz • • • , = jZ. In dem Erzeugendensystem hegt kein Vektor im Erzeugnis 
der anderen. Das Erzeugendensystem ist also „groß" und minimal, aber das daraus 
erzeugte Gitter besitzt das im Vergleich zu den pi relativ kleine Minimum j . 

Bei diesem Beispiel handelt es sich um einen atypischen Extremfall. Viele Git- 
ter, die in der Theorie der quadratischen Formen behandelt werden, etwa ganzzahlige 
Gitter, Gitter mit hoher Packungsdichte oder durch algebraische Codes konstruierte 
Gitter, besitzen eher große Minima. In disem Fall wird der Überdeckungssatz 3.14 
praktisch relevant und besagt, daß minimale Erzeugendensysteme solcher Gitter rela- 
tiv klein sind. 

Im folgenden untersuchen wir die praktische Anwendbarkeit von Satz 3. 14 und Al- 
gorithmus 3. Als Kern von Algorithmus 3 benutzen wir in unseren Versuchen den 
MLLL von Pohst für ganzzahlige Gitter. Bei nicht ganzzahügen Gittern verkom- 
pliziert sich die Situation lediglich bei der Vorbereitung des Inputs für unseren Algo- 
rithmus. Ist das Gitter zwar nicht ganzzahlig, aber rational, so können wir es durch 
geeignetes Skalieren auf ein ganzzahliges Gitter transformieren. Handelt es sich um 
ein echt reelles Gitter, so werden wir zunächst eine Approximation auf ein rationales 
Gitter vornehmen. Dieses Verfahren wird beispielsweise in [BP89] benutzt und seine 
Auswirkungen auf den LLL analysiert. Um die Problematik der mit einer Approxima- 
tion verbundenen Rundungsfehler hier zu vermeiden, betrachten wir nur ganzzahlige 
Gitter. 

Lenstra, Lenstra und Loväsz beweisen in [LLL82] für die Laufzeit des LLLs eine 
asymptotische obere Schranke von 0(n'^ log B) arithmetischen Operationen. Exper- 
imente zeigen, daß diese Schranke pessimistisch ist; dies gilt ebenfalls für die oben 
erwähnten Modifikationen des LLL wie den MLLL. Bei sehr vielen Inputvektoren 
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wird jedoch auch der MLLL merkbar langsam. Die naive Idee, sukzessive den ML- 
LL auf n + 1 Vektoren anzuwenden, führt zu keiner Verbesserung. Wir nennen diese 
Variante den inkrementellen MLLL. Die Abbildungen 3.3 und 3.4 zeigen die Laufzeit- 
en vom MLLL, dem inkrementellen MLLL und Algorithmus 3 im Vergleich. Wir 
wenden die Algorithmen auf Erzeugendensysteme mit Vektoren an, deren Koeffizien- 
ten ganzzahhg, pseudo-zufällig und vom Betrag kleiner als eine Konstante K sind. 
Insbesondere ist die Norm der Vektoren durch B = \fnK beschränkt. 




1 I 1 1 1 ^ ^ ^ 1 1 1 1 

100 200 300 400 500 600 700 800 900 1000 

number of generators 

MLLL, -I- incremental MLLL, □ Algorithm 3 

Abbildung 3.3: Laufzeiten in Dimension 20 mit K = 100 



Die Auswertung der Messungen in den Abbildungen 3.3 und 3.4 zeigt, daß die An- 
zahl der Updateoperationen in Algorithmus 3 sehr klein in Beziehung zur Größe des 
Erzeugendensystems ist. Typische Determinanten der Anfangsgitter von vollem Rang 
sind in der Größenordnung von 2^°°, die Determinanten der erzeugten Gitter sind 
sehr klein. Die Determinanten der Anfangsgitter enthalten oft sehr große Primteil- 
er, so daß die Anzahl der Updateoperationen sehr viel kleiner ist als die in Korol- 
lar 3.15 angegebene Schranke, bei der die vorkommenden Primteiler nach unten auf 2 
abgeschätzt werden. 

Die Abbildungen 3.3 und 3.4 stammen von Frank Vallentin, der auch das Testpro- 
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1000 



dimension = 20, ränge = [-1000, 1000] 




10 

100 200 300 400 500 600 700 800 900 1000 
number of generators 

MLLL, + incremental MLLL, □ Algorithm 3 

Abbildung 3.4: Laufzeiten in Dimension 20 mit K = 1000 

gramm entwickelt hat. Dabei wurde als Kern Victor Shoups NTL^ Bibliothek benutzt, 
aus welcher die rationale Arithmetik und der MLLL stammen. 



3.7 Implementierungshinweise 

Die oben beschriebenen Verfahren führen, wie an den Beispielen gezeigt, zu in der 
Praxis effizienten Implementierungen mit beweisbaren Komplexitätsabschätzungen. 
Dennoch sind mit einfachen Heuristiken und Methoden in vielen Fällen weitere erhe- 
bliche Laufzeitgewinne zu erzielen. Die Liste der folgenden Hinweise sollte nur als 
Wegweiser zu einer Implementation von schnellen und robusten Bibhotheksfunktion- 
en dienen. 



3.7.1 Das Zerlegungsproblem 

Beim Zerlegungsproblem ist eine Basis oder Grammatrix des Gitters bereits bekan- 
nt. Daraus kann man unmittelbar die Invarianten Rang und Determinante des Gitters 

*URL: http://www.shoup.net/ntl/ 
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sehr leicht berechnen. Diese können als Abbruchkriterium im Schleifendurchlauf von 
Algorithmus 2 verwendet werden. 

i) Wenn im Anschluß an Zeile 15 in Algorithmus 2 das Gitter Li denselben Rang 
wie das Gitter L besitzt, dann kann es nur ein triviales orthogonales Komple- 
ment besitzen, also ist bereits auch das Gitter L unzerlegbar und der Algorith- 
mus kann beendet werden. 

ii) Wenn im Anschluß an Zeile 15 in Algorithmus 2 die Determinante des Gitters 
Li0. . .®Lk gleich der Determinate von L ist, dann gilt Li®. . .®Lk = L, und 
wir haben eine Zerlegung von L gefunden, wenn auch noch nicht alle Vektoren 
in S abgearbeitet sind. 

Ist ein Gitter zerlegbar, so kann diese Eigenschaft gelegentlich schon an einer 
LLL-reduzierten Grammatrix des Gitters abgelesen werden. Wenn beispielsweise alle 
sukzessiven Minima nahe beieinander liegen, dann sind zerlegbare Vektoren relativ 
lang im Vergleich zu unzerlegbaren Vektoren. Der LLL-Algorithmus wird dann mit 
großer Wahrscheinlichkeit bereits eine Basis aus unzerlegbaren Vektoren finden. Die 
zu dieser Basis gehörige Grammatrix besitzt dann (bis auf eine Permutation der Ba- 
sisvektoren) die Gestalt einer Blockmatrix. 

3.7.2 Die Basisberechnung aus einem großen 
Erzeugendensystem 

Der Algorithmus 3 ist leicht parallelisierbar. Eine Erzeugendenmenge kann in kleinere, 
etwa gleich große Teilmengen zerlegt werden. Die Basen der aus den Teilmengen 
erzeugten Gitter können jeweils mit dem Algorithmus 3 nach dem „divide and conquer"- 
Verfahren parallel berechnet werden. In einem letzten Schritt werden die Basisvek- 
toren der Teilgitter wieder zu einem Erzeugendensystem S zusammengefügt, welches 
wiederum als Input für Algorithmus 3 dient. Diese Vorgehensweise empfiehlt sich 
allerdings nur für sehr große Erzeugendensysteme. In der Regel ist selbst eine häufige 
Ausführung der Lookup-Funktion in Zeile 8 schneller als der zusätzüche Aufwand 
durch wenige neue MLLLs. 
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Quadratische Formen und ihre Gitter finden in einem Gebiet der Informationstheorie 
wichtige Anwendungen. Bei der Übertragung von Informationen über ein störungsan- 
fälliges Medium ist es von eminenter Bedeutung, auf der Empfängerseite die gesende- 
ten Informationen korrekt und vollständig wieder herstellen zu können. Die grund- 
sätzhche Idee zur Fehlerkorrektur besteht in der Erhöhung der Redundanz der In- 
formationen, um kleinere Fehlübertragungen korrigieren zu können. Eine solche 
Möglichkeit ist die Verwendung fehlerkorrigierender Codes. Viele dieser Codes mit 
guten Korrektureigenschaften können aus Gittern konstruiert werden. Diese Meth- 
oden greifen aber erst, wenn die Informationen bereits als codierte, digitale Daten 
vorüegen. Bei der Übertragung durch Kabel oder Funkwellen erreichen die Informa- 
tionen den Empfänger zunächst analog und in der Regel gestört. Ein Analog-Digital- 
Wandler setzt die analogen Signale in digitale um. Schon in diesem Schritt werden 
kleinere Störungen korrigiert. Hierfür werden die gesendeten Daten nur aus einem 
kleinen, fest definierten Zeichenvorrat, dem sogenannten Code Book, gebildet. Ein 
empfangenes, gestörtes Datum wird mittels eines Quantizers, manchmal auch Ap- 
proximation oder Compressor genannt, auf das „nächstliegende" Zeichen im Code 
Book abgebildet. Bereits diese informale Beschreibung läßt die Beziehung zwischen 
Quantizem und Voronoiizellen erahnen. Eine Methode, Quantizer mit guten Korrek- 
tureigenschaften zu konstruieren, basiert auf Gittern im eukhdischen Raum. 

Definition 4.1. Sei P eine diskrete Punktmenge in E. Eine Funktion q : E ^ P heißt 
Quantizer, wenn für jedes v G E gilt, daß \ \v — q{v) \\ ^ | |x — p| | für alle p E P. 

Ein Quantizer ordnet also jedem Punkt im euklidischen Raum den nächstliegenden 
Punkt oder, falls dieser nicht eindeutig ist, einen der nächstliegenden Punkte in P zu. 

Ein Quantizer ist also im wesenthchen bereits durch die Angabe seiner Wertemenge 
charakterisiert. Lediglich für Punkte in E, die sich genau gleich weit von zwei Punk- 
ten in P befinden, ist der Funktionswert für zwei Quantizer von P nicht notwendiger- 
weise gleich. 

Man kann sich die Funktion eines Quantizers als hochdimensionalen Rundungsprozeß 
oder auch als eine verlustbehaftete Datenkompression vorstellen. Für die Qualität 
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eines Quantizers können je nach Anwendung verschiedene Maße dienen. Ein uni- 
verselles und weit verbreitetes ist der durchschnittliche quadratische Fehler je Dimen- 
sion 

Gg = - [ \\v-q{v)\\^ ■p{v)dv, (4.1) 
Je 

wobei die Funktion p die Wahrscheinschlichkeitsdichte der Datenquelle beschreibt 
([CS92], S. 58). 

Wenn man für P die Menge der Punkte eines Gitters L wählt, so läßt sich die 

Konstruktion von gitterbasierten Quantizern elegant mit dem geometrischen Konzept 
der Voronoizelle eines Gitters beschreiben. Im folgenden Abschnitt werden einige 
wichtige Definitionen und Eigenschaften der Voronoizellen von Gittern zusammenge- 
faßt. 

Alle Gitter L in diesem Kapitel haben vollen Rang n, d.h., sie spannen E als R- 
Vektorraum auf. 

4.1 Die Voronoizelle eines Gitters 

Zunächst wird der zentrale Begriff des Voronoibereichs eines Punktes eingeführt. 
Definition 4.2. Sei P eine diskrete Teilmenge von E. Für ein v E P heißt die Menge 

Vy{P) = {x & E \ \\x — v\\ ^ — w\\für alle w G P} 

der Voronoibereich von v in P. 

Ist V £ L ein Gitterpunkt, so nennt man Vu{L) die Voronoizelle von v in L. Wenn 
keine Verwechselung möglich ist, dann schreibt man auch kurz . 

Offensichtlich gilt q^^{v) C 14 (P) für ein v £ P. Im allgemeinen gilt jedoch 
keine Gleichheit, weil es Punkte x G E gibt, die sich in mehr als einer Voronoizelle 
befinden. 

Die Punkte in E, die sich am weitesten von Gitterpunkten entfernt befinden, kennze- 
ichnen eine wichtige Invariante des Gitters, den Überdeckungsradius. 

Definition 4.3. Für eine Teilmenge P C. E heißt die Zahl 

R{P) = sup inf IIa; - vi I G U {oo} 

der Überdeckungsradius von P. 

R(L) ist der größtmögliche Abstand, den ein Punkt in E von einem Gitterpunkt 
haben kann. Das Auffinden von Punkten mit dieser Eigenschaft, den sogenannten 



49 



4 Der Quantizer der Vorono'izelle eines Gitters 

tiefen Löchern ist hilfreich zur Konstruktion von Gittern mit hoher Packungsdichte 
(s. [CS92] Chapt. 6). 

In diesem Kapitel werden einige geometrischen Eigenschaften von Voronoizellen 
vorgestellt. 

4.2 Die Voronoizelle als Polytop 

Das erste Ziel ist es, die Voronoizelle Vq des Ursprungs eines Gitters als Polytop zu 
beschreiben und ihre Facetten zu bestimmen. 

In diesem Kapitel tritt der Aspekt der definierenden quadratischen Form in den 
Hintergrund. Zur Vereinfachung wird hier lediglich das euklidische Skalarprodukt 
{v,v) benutzt, das wegen der zahlreichen Klammern in einigen Formeln auch a\&v - v 
geschrieben wird. Falls v ^ E\0, schreiben wir H{v) = {x ^ E\x ■ v = v ■ v} für 
die Hyperebene mit dem Normalenvektor v und = {x ^ E \ x ■ v ^ v ■ v} für den 
dadurch bestimmten abgeschlossenen Halbraum, der die enthält. 

Wenn P eine diskrete Menge mit G P ist, dann kann man in dieser Notation den 
Voronoil^ereich von als Durchschnitt von Halbräumen beschreiben. 

Satz 4.4. Sei P eine diskrete Menge und G P. Dann ist 

veP\{Ö) 

Beweis. Sei 7^ G P und x e E. Dann gilt 

llx — 0|| ^ ||x — ■ü|| 

x-x^{x — v)-{x — v) = x- x — 2x-v + v- v 

x ■ V ^ ■ V 

X ■ ^ ■ 

X e m^. 

Damit folgt 

□ 

Im Fall, daß P ein Gitter ist, folgt aus der Definition 4.2 und der Eigenschaft P = 
P — vfür alle v G P unmittelbar das 

Lemma 4.5. Sei L ein Gitter auf E. 
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i) Für V E L gilt Vv = v + Vq. 

ii) Für a G 0{L) gilt Vbo" = Vq, inbesondere ist — Vq = Vq. 

Die Voronoizelle eines beliebigen Gitterpunktes ist ein Translat von Vq. Deshalb 
schreiben wir auch „die" Voronoizelle von L, wenn wir nicht eine ausgezeichnete 
Zelle meinen. 

Lemma 4.6. Es seien v,w E Vw, x eVw und v E L \ {w} mit \ \x — w\ \ = | |a; — ^1 1, 
dann ist auch x eVv. 

Beweis. Sei z E L \ {v,w}. Aus x E Vu, folgt, daß ||z — x|| ^ ||x — w|| = ||x — ?;||. 
Also ist X eVy. □ 

Bemerkung 4.7. Aus diesem Lemma folgt u.a. der wichtige Satz, daß T = {Vy \ v G 
L} eine Seite-an-Seite Pflasterung des Raumes E ist. 

Lemma 4.8. Sei C C. E konvex und abgeschlossen. Die Menge C ist genau dann 
beschränkt, wenn sie keinen Strahl enthält. 

Beweis. Wenn C einen Strahl enthält, dann ist C offenbar nicht beschränkt. Sei C 
nun nicht beschränkt und ohne Beschränkung der Allgemeinheit sei die Sphäre mit 
Radius 1 um den Nullpunkt 50(1) = {x e E \ \\x\ \ = 1} C C. Da die Menge C nicht 
beschränkt ist, enthält sie eine Folge {xn)n ^ C mit > n. Wegen der Kom- 

paktheit von 5*0(1) enthält die Folge (|]f^)n ^ 'S'o(l) nach dem Satz von Bolzano- 
Weierstraß eine konvergente Teilfolge ij^^^)nk ) deren Limes wir mit x bezeichnen. 
Wir zeigen nun, daß C den Strahl Ä+x enthält. Sei a G R'^. Für fast alle k ist Xn^. > 
a und wegen der Konvexität von C dann auch ^^ Xr^. G [0, Xn^.] C C fl 50(0;). 
Aus der Abgeschlossenheit von C folgt, daß limfe p^a^nfc = ax e C und damit die 
Behauptung. □ 

Lemma 4.9. Vq{L) ist beschränkt. 

Beweis. Seien w,x e E rmtw+ax G Vq für alle a G und Sx = {v G L \ x-v ^ 
0}. Für alle a G M^*' und v E Sx gilt, daß w + ax G Vq 44> {w + ax) ■ ^ ■ ^v. 
Damit folgt x • f = 0. L ist ein Gitter auf E und somit ist dim(5':c)]K = dimE* und 
folglich X = 0. Vb ist konvex und enthält keinen Strahl, womit ist gezeigt, daß Vq 
beschränkt ist. □ 
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4.3 Voronoirelevante Vektoren 

Definition 4.10. Sei v e L\0. 

i) Wir nennen v einen Hüllenvektor, wenn H{lv)nVo ^ ^,also seine Hyperebene 
die Vorono'izelle berührt. 

ii) V heißt voronoirelevant oder kurz relevant, wenn H{-^v) einen n — 1-dimensio- 
nalen (affinen) Schnitt mit Vq besitzt, d. h., H{^v) eine Wand von Vq enthält 
[ CS92 ]. Ist V relevant, so schreiben wir Fyßrdie Wand H{\v) n Vq. 

Die Menge der voronoifrelevanten Vektoren genügt, um Vq zu konstruieren. 

Satz 4.11. Die Wände Fy von Vq sind zentralsymmetrisch mit Mittelpunkt |v. 

Beweis. Sei t E Fy. Wie in 4.5 gezeigt, liegt —t + vinVy und wegen {—t + v)- = 
■ auch in H{^v), also mit Lemma 4.6 in Vq H H{^v) = Fy. Damit ist auch 
die Konvexkombination + ^{—t + v) = in Fy und Fy ist zentralsymmetrisch 
bezüglich ^v. □ 

Wir wissen nun, daß Vq als Polytop ein endlicher Durchschnitt von Halbräumen ist. 
Bereits 1908 hat Voronoi eine Kennzeichnung dieser Halbräume angegeben: 

Satz 4.12 ([Vor08]). Ein Vektor ^ v £ L ist genau dann relevant, wenn ibv die 
einzigen kürzesten Vektoren in der Nebenklasse v + L/2L sind. 

Beweis. Wegen der Wichtigkeit dieses Satzes und der Kürze des Argumentes geben 
wir hier den Beweis aus [CS92], S. 475^ an. 

Es seien v^w G L mit v — w E 2L,v ^ ±w und w ■ w ^ v ■ v. Dann sind 
t = ^{v + w) und u = ^{v — w) beide in L \ und man verifiziert x ■ v = x ■ {t + 
u) ^ ■ t + i( • ^ \{v ■ V + w ■ w) ^ \v ■ V. Damit ist v nicht relevant, weil 

X e HtCiHu^ X e Hy. 

Sei nun v G L nicht relevant und einer der kürzesten Vektoren in der Nebenklasse 
V + 2L. Dann gibt es ein w G L \ mit w ^ v und ■ ^ ■ ^w. Wir erhalten 

{v — 2w) ■ {v — 2w) = V ■ V — 4:V ■ w + Aw ■ w ^ V ■ V 

und haben mit v — 2w wegen ^ v — 2w ^ ±v einen weiteren kürzesten Vektor in 

+ 2L gefunden. □ 

'in der ersten Auflage fehlt wegen eines Druckfehlers die Formulierung dieses Satzes. In der zweiten 
Auflage ist der Satz enthalten, aber der Korrekturhinweis im Vorwort, S. xxviii, verweist auf die 
falsche Seite. 
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An dieser Stelle sei nochmals auf das Korollar 3.6 auf Seite 33 verwiesen, das genau 
die spitz unzerlegbaren Vektoren als voronoifrelevante identifiziert. 

Korollar 4.13. Für ein Gitter L gilt: 

i) Es gibt höchstens 2 • (2" — 1) relevante Vektoren. 

ii) Vo{L) ist endlicher Durchschnitt von Halbräumen, also ein Polytop, insbeson- 
dere konvex. 

iii) R{L) < oo. 

Bemerkung 4.14. Bereits Voronoi zeigte, daß die Schranke von 2 • (2" — 1) relevanten 
Vektoren scharf ist [Vor08, Vor09]. 

4.4 Der Quantizer eines Gitters 

Im Falle einer Datenquelle mit uniform verteiltem Output und einem Gitterquantizer 
vereinfacht sich die Gleichung (4.1). Dies motiviert die folgende 

Definition 4.15 ([CS92], Chapt. 2.3, Chapt. 21.). Für ein Gitter L von vollem Rang, 
heißt 



das normalisierte dimensionslose zweite Trägheitsmoment von L. 

Mathematisch exakt läßt sich diese Invariante von L im allgemeinen nur schwer 
errechnen. In [CS92], Chapt. 21 findet sich die Herleitung für einige wichtige Stan- 
dardpolytope, unter ihnen auch das Simplex. Letzteres wird dort benutzt, um mit 
Hilfe der Operation der Weyl-Gruppe den Quantizer G{L) für alle Wurzelgitter L zu 
berechnen. 

Für Gitter mit weniger Strukturinformationen ist es praktikabel, das Trägheitsmo- 
ment in Gleichung (4.2) durch numerische Monte-Carlo-Integration zu approximieren. 
Es ist jedoch nicht unmittelbar möglich, eine uniforme Verteilung von Samples in Vb 
zu erzeugen, weil die genaue Struktur von Vq im allgemeinen nicht bekannt ist. Man 
kann sich hier die Translationseigenschaft von Voronoizellen zunutze machen. Man 
definiert eine Abbildung Vq [0, 1)" folgendermaßen. Für eine feste Basis B von 
L bestimme man für jedes y € Vq den Koordinatenvektor u = bezüglich B 
und bestimme komponentenweise den nicht ganzzahligen Anteil x = {xi, - ■ ■ , Xn) = 
{ui — [uij , ■ ■ ■ ,Un— lun\)- Mit Hilfe der Translationseigenschaft von Voronoizellen 




(4.2) 
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verifiziert man leiclit, daß diese Abbildung bijektiv ist. Somit kann man mit ein- 
er uniformen Datenquelle in [0, 1)" eine uniforme Verteilung in Vq erzeugen.^ Die 
Transformation mit B muß im Integral in Gleichung (4.2) durch eine entsprechende 
Substitution berücksichtigt werden. 

Im folgenden Algorithmus 4 bezeichnet S das arithmetische Mittel der Abstände 
der Samples zu Gitterpunkten, G{L) das approximierte Trägheitsmoment, = 
S)^ die Streuung und a = den Schätzwert des Fehlers. Für den korrekten Rück- 
gabewert müssen & und G{L) jeweils mit -(det L)"« skaliert werden. 



1: Input: Ein Gitter L Q E, gegeben durch eine Basis B und t EN, die Anzahl der 

zu berechnenden Samples. 

2: Output: G{L) und ä, der Schätzwert des Fehlers von G{L). 

3: // 1. InitialisieruDg. 

4: Initialisiere ein Array S'[O...A;] für die zu berechnenden Samples. 

5: // 2. Berechnung der Samples. 

6: for A; = to t do 
7: Wähle Xk G [0, l)*^ randomisiert. 
8: yk closest_vector(L, B ■ Xk) 
9: S[k]^\\yk-B-Xk\\'' 

10: end for 

11: // 3. Ergebnisaufbereitung. 

12: S 'S- i Yä //Arithmetisches Mittel über die S[k] 

15: G{L) ^ i(detL)-^ -S 
16: retum {G{L),a) 

Algorithmus 4: Numerische Approximation von G{L) 

Algorithmus 4 kann leicht verallgemeinert werden, um den mittleren quadratischen 
Fehler für Quantizer von beUebigen Vereinigungen von Nebenklassen eines Gitters 
Lc UcgC L + c für C = {ci, - ■ ■ ,Cm} C E zu berechnen. Hierzu bestimmt man 
nicht den Abstand von B ■ Xk zu L, sondern jeweils den zu L + c, c G C und setzt 
S[k] auf das Minimum aller berechneten Abstände. Für die Bestimmung von G{Lc) 
und des entsprechenden Schätzwertes des Fehler muß dem Umstand Rechnung ge- 
tragen werden, daß sich die Dichte der Gitterpunkte um den Faktor m gegenüber der 
Dichte von L vervielfacht. Dies wird durch Korrekturfaktor bei der Determinante 

^Die zugnmdeliegende geometrische Idee wird in [CS82b, CS82a] beschrieben. 



54 



4 Der Quantizer der Voronoizelle eines Gitters 

berücksichtigt. Dem liegt die Identität (^cgc* vol Vc)^ = dctL zugrunde, insbeson- 
dere ist es nicht notwendig, die exakten Volumina , • • • , zu kennen. Somit 
erhält man Algorithmus 5. 

1: Input: Ein Gitter L C E, gegeben durch eine Basis B, eine Menge C = 

{ci, • ■ ■ ! Cm} C E und t G N, die Anzahl der zu berechnenden Samples. 

2: Output: G{L) und a, der Schätzwert des Fehlers von G{L) 

3: // 1. Initialisierung. 

4: Initialisiere ein Array 5[0...fc] für die zu berechnenden Samples. 

5: // 2. Berechnung der Samples. 

6: for A; = to t do 

7: Wähle Xk G [0, 1)" randomisiert. 

8: Uk <— closest_vector(L, B ■ Xk) 

9: S[k] ^ Wvk - B ■ Xkll"" 

10: for j = 1 to m do 

11: Uk closest_vector(L, B ■ x^ — cj) 

12: S[k] ^ mm{S[k] ,\\yk - B ■ Xk + CjW^) 

13: end for 

14: end for 

15: //3. Ergebnisaufbereitung. 

16: S ^ i ^[i]. // Arithmetisches Mittel über die S[k] 

18: a^i(detL/|q2)-i.^ 

19: G(Le)^i(detL/|Cn-^.S 
20: retum {G{Lc), ä) 

Algorithmus 5: Numerische Approximation von G{Lc) 



4.5 Gitter mit kleinem Quantizer 

Das dimensionslose zweite normalisierte Trägheitsmoment eines Quantizers G{L) ist 
eine Invariante des Gitters L, die bis vor kurzem wenig untersucht wurde. Aus der 
Definition von G{L) in Gleichung (4.2) folgt unmittelbar, daß für ein festes Volu- 
men von Vb (L) das Trägheitsmoment minimal ist, wenn Jy^x ■ x dx minimal für alle 
Gitter L ist. Nach dem Isoperimeterprinzip besäßen n-dimensionale Kugeln das klein- 
ste Trägheitsmoment. Kugeln sind jedoch nicht raumfüllend wie Voronoizellen. Zur 
Zeit hat man kein anschauliches Verständnis von den Parametern, die ein niedriges 
Trägheitsmoment bestimmen. In einer neueren Arbeit erzielten Agrell und Eriksson 



55 



4 Der Quantizer der Voronoizelle eines Gitters 



[AE98] eine Reihe von zum Teil kontraintuitiven Ergebnissen über Gitterquantizer. 
Sie adaptieren dort einen Trainingsalgorithmus für Quantizer auf die Gittersituation. 
Die mit dieser Methode gefundenen Gitter befinden sich bereits relativ nahe an der 
von Conway und Sloane vermuteten theoretischen unteren Schranke für den mit- 
tleren quadratischen Fehler eines Gitterquantizers [CS 85], insbesondere verbessern 
sie die bekannten Rekorde in den Dimensionen 9 und 10, s. [CS92], S. 61. In Dimen- 
sion 10 fanden sie in Df^ ein Gitter, dessen sehr niedriges Trägheitsmoment sie mit 
0.070814 ± 0.000001 angeben. Der Umstand, daß ein relativ bekanntes Gitter wie 
D^Q in dieser Hinsicht bisher immer übersehen wurde, legt nahe, daß Gitterquantizer 
bis dahin nur oberflächlich untersucht wurden, ganz im Gegensatz zu den klassis- 
chen Problemen wie dem Packungs- und dem Überdeckungsproblem. In Dimension 
9 berechneten sie ein nichtklassisches Gitter mit der folgenden Erzeugermatrix 
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0.573/ 



und einem Trägheitsmoment von 0.071626 ± 0.0000002 (nach [AE98]), welches sich 
schon relativ nahe am von Conway und Sloane vermuteten theoretisch möglichen 
Minimum 0.070902... befindet. Die ersten acht Zeilen erzeugen das Gitter Dg. Die 
ersten acht Koordinaten des neunten Vektors beschreiben ein „tiefes Loch"^ von Dg. 
Fügt man das deep hole als Erzeuger zu Dg, hinzu, so erhält man das Gitter Eg, . Die ne- 
unte Koordinate mit dem seltsam anmutenden, numerisch approximierten Wert 0.573 
läßt sich geometrisch nicht einordnen. Das Minimum von Dg ist gleich 2, das Gitter 
von Agrell und Eriksson besitzt Minimum (2 • 0.573)^ = 1.313316. Die Voronoizelle 
dieses Gitters ist also relativ flach, während man nach der obigen Motivation eher eine 
kugelförmigere Voronoizelle erwartet hätte. 

Mit tn kann man in verschiedenen Dimensionen eine Reihe von Gittern finden 
mit fast so niedrigen Trägheitsmomenten wie die niedrigsten bekanntesten. Auch in 
höheren Dimensionen, als sie von Agrell und Eriksson veröffentlicht wurden, findet 
man mit tn Gitter mit Trägheitsmomenten in der Nähe der vermuteten theoretischen 
unteren Schranke, wie zum Beispiel G{A^^) ^ 0.070422 ± 0.000012. 



'Siehe S. 50. 
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